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Introduction 



La mecanique des fluides est une discipline ancienne, d'applications tres variees et 
encore en pleine evolution. II convient de toujours garder a l'esprit que revolution 
de cette discipline a eu deux moteurs, fortement imbriques : 

- l'explication des phenomenes naturels : les vagues, le vent, la force de resistance 
sur un corps en mouvement dans Fair ou l'eau, l'aspiration d'une cheminee, le 
mouvement des bulles, la chute d'objets leger (les feuilles des arbres), les vibrations 
provoquees par un ecoulement. . . 

- l'exploitation des fluides a des fins pratiques : fabrication d'embarcations, pompage 
de puits, adduction d'eau, application "energetiques" (moulins a eau ou a vent) 

La mecanique des fluides a cet avantage qu'elle fait partie de notre quotidien, et il 
est toujours bon d'apprehender un ecoulement de fluide tout d'abord avec sa seule 
intuition. Les equations de la mecanique des fluides ont une structure mathematique 
complexe, et doivent etre vues comme un ultime recours pour decrire ou quantifier 
un phenomene, la ou l'intuition s'arrete. Elles ne sont pas la mecanique des fluides, 
elles la decrivent. 

A titre de preliminaire a ce cours, on pourra par exemple se poser les questions 
suivantes : 

- pourquoi les portes claquent dans un courant d'air ? 

- pourquoi les bulles remontent dans l'eau ? 

- pourquoi un avion ne tombe-t'il pas sous son propre poids? 

- comment faire monter un fluide d'un point bas a un point haut ? 

- la quantite d'eau sortant d'un tuyau par unite de temps est-elle egale a celle qui 
y rentre ? 

- l'eau sortant d'un robinet a-t-elle la meme allure selon la valeur du debit ? 

- que faut-il faire pour siphonner un recipient ? 

- comment fonctionne une centrale hydroelectrique ? 

Toutes ces questions sont generalement posees plutot par les enfants, parce qu'ils 
cherchent a comprendre le monde qui les entoure. Avec l'age, la tendance intellec- 
tuelle naturelle est plutot a l'acceptation du phenomene. II est done important de 
questionner a nouveau sa perception des phenomenes physiques avant de vouloir 
les modeliser. Completez done cette liste de questions par les votres en continuant a 
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etudier ce cours. Le cadre forcement restreint de ce dernier ne permettra pas malheu- 
reusement de repondre a toutes, mais cette demarche intellectuelle est importante. 
Pensez egalement que les scientifiques qui ont etabli les equations de la mecanique 
des fluides sont alles de l'intuitif vers le modele mathematique, et non pas l'inverse. 

Quelques conseils de lecture de ce polycopie : 

- le chapitre 1 decrit ce qu'est un fluide et comment on peut le definir par un milieu 
continu. Sa lecture, bien que conseillee, n'est pas indispensable pour comprendre 
les chapitres suivants. 

- le chapitre 2 decrit ce qu'est une equation de bilan et definit mathematiquement 
le terme de flux convectif. On pourra eventuellement survoler ce chapitre a l'usage 



du debutant, et ne retenir que son resultat essentiel (2.6 ). On y definit egalement ce 
volume particulier qu'est le tube de courant. On y expose, ensuite, les equations de 
conservation pour un fluide sans detainer les forces exercees, et les trois equations 
de bilan fondamentales pour un fluide sont retrouvees : la masse, la quantite de 
mouvement et l'energie. Ce n'est finalement qu'un rappel des grands principes de 
la physique. Sa lecture est indispensable. 

- Le chapitre 3 peut etre lu indifferemment avant ou apres les trois precedents. II 
decrit les forces de volume et de contact exercees sur un fluide. On y deduit les 
lois de l'hydrostatique. La notion de force de frottement visqueuse y est presentee 
et le nombre de Reynolds introduit. Sa lecture est indispensable. 

- Le chapitre 4 est le coeur de ce document : sur la base des chapitres 2 et 3, les 
equations du mouvement d'un fluide sont presentees, sous forme globale, dans le 
cas particulier d'une machine fluide, et enfin sous forme locale. Aucune hypothese 
restrictive n'y est effectuee de telle sorte que les equations presentees sont appli- 
cables a tout type de fluide ou d'ecoulement. L 'hypothese incompressible y est 
discutee. Lecture indispensable. 

- Le chapitre 5 presente le modele du fluide parfait uniquement dans le cas incom- 
pressible. La formule de Bernoulli y est presentee, ainsi que la notion de charge. 
Les limites de cette formule et du modele de fluide parfait sont egalement presen- 
tees, ainsi que les differents cas de gains ou pertes de charge dans un ecoulement 
en tuyauterie. Lecture indispensable. 

- Le chapitre 6 introduit le caractere visqueux des fluides. Le modele de fluide dit 
"newtonien" est introduit sur la base de resultats experimentaux sur l'experience 
de Couette. Les equations de Navier-Stokes en sont deduites et adimensionalisees 
pour faire apparaitre des nombres adimensionnels, en particulier le nombre de Rey- 
nolds. On detaille ensuite l'ecoulement de Poiseuille, avec les importantes conse- 
quences sur la perte de charge en ligne et les modeles de milieux poreux. Lecture 
indispensable. 

- Le chapitre 7 (en cours d'ecriture) traite le cas particulier des ecoulements dits 
"rampants", c'est-a-dire a tres faible nombre de Reynolds, avec pour application 
le calcul de la trainee sur une sphere. 

Ce polycopie est inspire de plusieurs presentations differentes de la mecanique des 



fluides. Citons tout d'abord le livre passionnant de Guyon et al. (2001) qui est a 



notre connaissance l'ouvrage le plus exhaustif sur le sujet (et en plus en frangais!). 
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II est notamment le seul a traiter en detail les fluides non-newtoniens et les effets 
dans des referentiels tournants. Attention, la seconde edition (couverture verte) est 
bien plus complete que la premiere (couverture rouge). Une excellente introduction 
a la turbulence est egalement proposee. 



On pourra consulter l'excellent livre de Chassaing (2000) (chercheur a l'lnstitut de 



Mecanique des Fluides de Toulouse), tres complet notamment dans le traitement 
des couches limites, et d'une rigueur mathematique remarquable . 



Le livre de Candel (1995), outre sa presentation irreprochable, traite en detail cer- 



taines applications inusuelles (eolienne, fusee,. . .). II partage avec l'ouvrage de White 



(1994) et ce polycopie une presentation preferentielle des equations sous forme volu- 



mique. Le livre de White (1994) presente par ailleurs un grand nombre de donnees 



technologiques, d'applications pratiques et d'exercices originaux. C'est une source 
d'information incontournable (mais en anglais) pour tout probleme pratique. 



Chapitre 1 

Description d'un fluide 



1.1 Qu'est-ce qu'un fluide ? 



En relation avec le cours de thermodynamique, nous incluons dans la notion de 
fluide les liquides et les gaz. On pourrait tout d'abord demander au lecteur, avant 
d'aller plus loin, de se questionner sur la definition d'un fluide : le sens commun 
suggere qu'un fluide, en opposition avec un solide qui sait maintenir sa forme tout 
seul, necessite un contenant, sans quoi il "coule" (pour un liquide) ou "s'echappe" 
(pour un gaz). 

Concernant les gaz, on sait qu'ils sont constitues d'atomes ou de molecules en mou- 
vement (dit d'agitation thermique) qui s'entrechoquent en permanence. Cette des- 
cription a ete proposee a l'origine par BoltzmamQ et conduit a la fameuse theorie 
cinetique des gaz. Si le gaz est suffisamment dilue (limite du gaz parfait), les mo- 
lecules n'interagissent pas autrement que par des chocs elastiques, car elles sont en 
moyenne suffisamment eloignees les unes des autres pour que les forces electrosta- 
tiques soit negligeables. 

Les liquides nous sont beaucoup plus familiers que les gaz [^] , mais leur structure est 
pourtant plus complexe. lis sont l'intermediaire entre le solide et le gaz. Le solide 

1 Ludwig Boltzmann (1844-1906) : il est le fondateur de la description microscopique des gaz 
a partir des principes de la mecanique classique. II fut notamment le premier a donner une in- 
terpretation microscopique de l'entropie et de la notion d'irreversibilite. Ses idees audacieuses, 
qui constituent le fondement de la thermodynamique moderne, furent vivement contestees de son 
temps, notamment par Mach et Ostwald, pour qui la science doit se baser sur des observations 
macroscopiques et ne pas faire appel a des entites microscopiques (l'atome) dont l'existence ne 
peut etre observee. Decourage, il ecrivit en 1898 : "je suis conscient de n'etre qu'un faible individu 
luttant contre le courant de l'epoquc". Par ailleurs de tendance depressive, il commit deux tenta- 
tives de suicide, dont la seconde lui fut fatalc, en 1906. II ne connut ainsi jamais les consequences 
de sa theorie. La preuve de l'existence de la structure atomique de la matiere fut apportee deux 
ans apres sa mort. 

2 Tout d'abord parce qu'on les "voit", et que, contrairement aux gaz, on ressent immediatement 
leur contact avec la peau (signification du mot "mouille"). L'air est pourtant en permanence en 
contact avec notre peau et lui applique une force de 1 kg par cm 2 ! 
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1.2 Proprietes. 



maintient ses atomes fixes, enfermes dans une "cage" constitute par les atonies voi- 
sins, a l'exception de vibrations de faible amplitude d'origine thermique. A l'inverse, 
le gaz contient des atomes ou molecules libres les uns par rapport aux autres, n'in- 
teragissant que par des chocs. On retrouve les deux aspects dans un liquides : les 
molecules sont tres peu mobiles, enfermees dans des cages formees par leur voisines, 
mais ces cages sont temporaires, et, en presence d'une force suffisante, le mouvement 
relatif des molecules moyen peut devenir non nul. C'est ainsi qu'un liquide s'ecoule. 

Pour certaines matieres, il est parfois difficile de faire la difference entre un compor- 
tement liquide et solide. Par exemple le sable peut s'ecouler alors qu'il est constitue 
de particules solides. Pire encore, si l'on tape fortement sur du sable mouille, on 
ressent un comportement solide, mais on peut enfoncer lentement sa main dans le 
sable comme dans un liquide. Les glaciers apparaissent solides a notre echelle, mais 
s'ecoulent sur une echelle de temps de l'ordre du siecle. Certaines pates (dont cette 
fameuse gelee translucide vendue pour amuser les enfants. . .) apparaissent comme 
solides lorsque la sollicitation est rapide, mais coulent si on l'etire lentement. Lors 
du pressage de feuilles metalliques a froid, le solide prend la forme de l'outil et ne 
revient p£lS 8b Set forme initiale (heureusement !) car on depasse sa limite d'elasticite : 
pendant un temps tres court, il coule (ou parle de comportement plastique). 

II est bon de connaitre au moins l'existence de ces fluides "bizarres", qui ne le sont 
en realite que parce que leur comportement depend de l'echelle de temps consideree. 
L'etude du comportement d'une substance sous une sollicitation donnee porte le 
nom de rheologie. Nous nous limiterons a priori dans ce cours a des cas "courants", 
bien que certains principes enonces soient d'une portee plus generale. 



1.2 Proprietes. 

1.3 Description comme un milieu continu. 
1.3.1 Separation des echelles 

Nous avons mentionne ci-dessus le caractere "granulaire" des fluides : ils sont consti- 
tues d'entites elementaires, atomes ou molecules. On congoit qu'il serait difficile de 
prendre en compte les mouvements individuels de toutes ces entites, pour etudier 
le mouvement d'un liquide dans un tuyau ou de l'air autour d'une aile d'avion. . . 
Ce serait de plus contre-intuitif puisqu'a l'echelle humaine, nous ne voyons ni ne 
sentons pas individuellement les molecules |^J mais nous savons en revanche appre- 
hender leur mouvement (ou leur non- mouvement) d'ensemble a notre echelle (dite 
macroscopique). 

3 I1 est interessant de noter du reste que la mecanique des fluide est bien anterieure a la descrip- 
tion atomistique de la matiere, et que les premieres equations descriptives furent de fait etablies a 
l'echelle macroscopique. 



1.3 Description comme un milieu continu. 
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L'idee de l'approche macroscopique est basee sur le fait que des molecules suffi- 
samment voisines ont le meme mouvement d'ensemble, et ne different que par des 
mouvements d'agitation thermique aleatoires^J Dans ces conditions, il est possible 
de trouver un petit volume spatial e, contenant suffisamment de molecules pour 
pouvoir effectuer une moyenne de grandeurs microscopiques sur ce volume, mais 
suffisamment petit pour ne pas gommer les variations macroscopiques. 

Prenons une analogie : imaginons un embouteillage s'etendant sur plusieurs dizaines 
de kilometres que nous regardons depuis un helicoptere R et examinons ce que nous 
voyons en fonction de la hauteur de l'helicoptere (Fig. |f .1 [ ) : 

- si nous sommes proches du sol, nous verrons les deplacements individuels de 
quelques voitures, mais sans voir revolution d'ensemble. Nous avons une vue 
microscopique (Fig. |l.l[ c) 

- a l'inverse si nous volons tres haut, nous verrons plusieurs milliers de voitures et les 
variations de densite de l'embouteillage (oil encore les variations de distances entre 
deux voitures successives) deviendront visibles : dans certaines zones les voitures 
seront pare-chocs contre pare-chocs et a 1 'arret, dans d'autres elles seront un plus 
ecartees les unes des autres et en mouvement. C'est une vue macroscopique 
(Fig. [TTTja). Notons que les variations de densite apparaissent alors comme une 
fonction continue, voire derivable^] 

- Si nous volons a une altitude intermediaire entre ces deux points de vue, nous 
voyons un grand nombre de voitures (six sur la figure) equidistantes et se deplagant 
en gros a la meme vitesse : c'est le volume e dont nous parlions. II est tres grand 
a l'echelle de la voiture, mais tres petit a l'echelle de l'embouteillage. On parle de 
volume mesoscopique (meso = du milieu, scopos = observer) ou de meso-echelle 

(Fig-igb). 

Si cette separation des echelles est possible^ nous pouvons assimiler e a une 
variation infinitesimale dx de longueur macroscopique x et definir la densite locale 



de voitures par une fonction continue p(x) (la courbe tracee sur la figure f.f ) telle 
que le nombre de voitures a l'interieur de dx est : 



dN = p(x) dx 

4 Avertissons le lecteur des maintenant que ces mouvements thermiques sont responsables des 
phcnomenes de transfert de chaleur et de matiere, ainsi que du frottement visqueux. L'approche 
macroscopique ne consiste done surtout pas a les negliger, mais plutot a les moyenner pour obtenir 
des lois sur les grandeurs macroscopiques. On retrouve ainsi les lois de Fourier, Fick, ainsi que la 
loi de frottement visqueux des fluides dit newtoniens. Mathematiqucment ce moyennage est tres 
complexe car les interactions deux a deux entre molecules doivent etre considerees. 

5 L'analogie du trafic routier avec les fluides n'est pas totalcmcnt fortuitc. II cxistc effectivement 
des modeles de milieux continus pour ce genre de probleme, semblables d'une certainc manicrcs 
aux equations du mouvement des fluides compressibles. 

6 Ce n'est plus vrai en toute rigueur dans le cas important des phcnomenes de chocs. 

7 Elle ne Test pas toujours en mecanique des fluides, ou bien dans le cas de milieux trop rarefies, 
ou encore dans les phcnomenes de chocs, ou les grandeurs macroscopiques varient sur quelques lon- 
gueurs moleculaires. Ce type d'hypothese est egalement utilise pour decrire les milieux diphasiques 
(liquides a bulle, suspensions,. . .) comme des milieux monophasiques. 
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1.3 Description comme un milieu continu. 



p(x) 



a. 



b. 



c. 



1 

^- 

e 
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Fig. 1.1 - Un embouteillage vu depuis trois echelles tres differentes : celle de la 
voiture (en bas), celle de l'embouteillage (en haut), et une echelle intermediaire 
contenant beaucoup de voitures, mais oil la densite varie peu. 

de telle sorte que le nombre de voitures entre les abscisses a et b s'ecrit : 



N — p(x) dx 



1.3.2 Definition de la masse volumique 



C'est ce type d'approche que l'on utilise pour decrire un fluide par un milieu continu 
(le volume intermediaire contient dans ce cas plus de 6 "voitures"), sauf que nous 
avons un milieu tridimensionnel. Le volume e est grand par rapport a la molecule (il 
en contient plein), mais petit par rapport a la longueur caracteristique du phenomene 
(Figure 1.2) : ce sera done un element infinitesimal de volume macroscopique dV = 



1.3 Description comme un milieu continu. 



dx dy dz centre sur le point (x, y, z). 
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Fig. 1.2 - De gauche a droite : echelle macroscopique, mesoscopique, microscopique. 

La masse de toutes les molecules mj contenues dans e est aussi, a l'echelle macro- 
scipique, une quantite infinitesimale wdM. On definira done la masse volumique 
locale par : 

rrii dM 

p(x,y,z,t) = ^ r ^- = — (1.1) 

dont l'unite S.I. est le kg/m 3 . Notons que p depend a priori du temps puisque des 
molecules peuvent entrer et sortir du volume dV. Si nous cherchons la masse M(t) 
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1.3 Description comme un milieu continu. 



d'un volume V macroscopique (celui de la figure 1.2 gauche, par exemple), il suffit 
d'integrer cette relation sur tout ce volume : 




Remarque : nous prendrons pour habitude de ne pas ecrire la dependance des in- 
tegrandes en (x,y,z,t), afin d'alleger les ecritures, mais il convient de retenir que 
ceux-ci dependent a priori des 3 variables d'espace, et ne pas les sortir par erreur 
des integrates . 

Un cas particulier important sur lequel nous reviendrons souvent est celui oil la masse 
volumique peut etre considered constante a la fois dans l'espace et le temps. On parle 
alors de fluide incompressible. Dans ce cas la masse d'un volume V est simplement 
donnee par M = pV. Dans la plupart des applications, les liquides seront considered 
incompressibles. Les gaz en ecoulement peuvent etre aussi tres souvent considered 
incompressibles. En fait, la vazlidite de cette approximation ne porte pas sur le type 



du fluide, mais sur la valeur du nombre de Mach (voir section 4.4.4). 



1.3.3 Vitesse et quantite de mouvement 

Comment maintenant definir la vitesse macroscopique du fluide en un point ? Notre 
but est d'etudier l'ecoulement des fluides sous l'influence de forces, qui modifient 
la quantite de mouvement P du fluide. II est done pertinent de definir la vitesse 
moyenne a partir d'une moyenne de la quantite de mouvement, soit : 

v (x, y, z, t) = -=k = -Tfr, (1-3) 

2^ rrii pdV 

oil Ton a utilise le fait que, d'apres la section precedente, la masse du volume de 
moyennage e est ^ = M = pdV 



La quantite de mouvement globale d'un volume V s'ecrit : 

P(*)= fffpvdV (1.4) 



v 



1.4 Grandeurs locales et globales 
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On notera que v», vitesse des molecules, inclut a la fois le mouvement macroscopique 
du fluide, mais aussi l'agitation thermique microscopique. L'operation de moyennage 
de la quantite de mouvement permet de faire disparaitre cette derniere contribution. 
Ainsi dans un fluide au repos, v = mais v, ^ 0. 



1.3.4 Grandeurs energetiques 

Par un processus de moyennage similaire, il nous est possible de moyenner l'ener- 
gie cinetique du fluide. La mecanique d'un ensemble de masse ponctuelles fournit le 
cadre theorique pour ce faire. Nous ne souhaitons pas ici entrer dans le detail, mais il 
nous parait important de rappeler un resultat essentiel sous-tendant le premier prin- 
cipe de la thermodynamique. L'energie cinetique dans le referentiel du laboratoire 
d'un ensemble de particules en mouvement est egale a : 

- l'energie cinetique de leur centre de masse, oil serait concentree la somme des 
masses des particules. Cette contribution correspond a l'energie cinetique macro- 
scopique du fluide en ecoulement notee K. Pour le volume dV de la figure |1.2[ on 
a done dK = ^dm v 2 = \pdVv 2 . L'energie cinetique macroscopique d'un volume 
de fluide V s'ecrit done : 

K(t) = JJJ \ P v 2 dV (1.5) 
v 

- l'energie cinetique des particules dans le referentiel du centre de masse. C'est 
l'energie d'agitation thermique du fluide, qui contribue a l'energie interne 

Pour cette derniere, il est plus pratique d'utiliser une grandeur massique. Ainsi, pour 
un gaz parfait (c'est aussi vrai pour un liquide incompressible, cf. cours de ther- 
modynamique), l'energie interne massique est simplement definie par du = C v dT. 
L'energie interne par unite de volume est obtenue en multipliant l'energie massique 
u (en J/kg) par la masse volumique p (en kg/m 3 ), et l'energie interne d'un volume 
de fluide V sera done donnee par : 



U(t) = pudV (1.6) 




1.4 Grandeurs locales et globales 



A ce point, nous avons defini le fluide comme un milieu continu, oil toutes les "as- 
perites" microscopiques sont gommees. Nous savons definir la masse, la quantite de 
mouvement et les differents types d'energies d'un volume de fluide V . 

8 En fait C'EST l'energie interne en Fabsence d 'interactions entre particules. Dans le cas inverse, 
l'energie interne comporte aussi une contribution de l'energie potentielle d'interaction. 
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1.5 Volume fixe ou mobile? 



Notons que toutes les grandeurs macroscopiques que nous avons definies ci-dessus 
(M, P, U et K) sont des grandeurs extensives, c'est-a-dire proportionnelles a la 



quantite de matiere. Ainsi, les relations (1.2), 
forme : 



(1.4), (1.6) et (1.5) sont toutes de la 



G(t) = JJJgdv, 

V 

oil g est une grandeur volumique, exprimee en (Unite de G) / (Unite de volume). 
Les grandeurs volumiques associees a M, P, U et K sont respectivement p, pv, pu 
et pv 2 /2. Le but des chapitres suivants est de calculer les variations temporelles de 
ces grandeurs sur la base de grands principes physiques. 

II convient de bien comprendre que G est une grandeur globale, calculee pour un 
volume macroscopique V donne, et ne depend que du temps. En revanche g est 
une grandeur volumique locale, dependant en plus du point macroscopique (x, y, z) 
considere. Mathematiquement G est une fonction d'une variable et admet une derivee 
temporelle, alors que g est une fonction de 4 variables et admet des derivees partielles 
par rapport a ces 4 variables. 



1.5 Volume fixe ou mobile? 



Notons que dans les definitions ci-dessus, aucune restiction sur le volume V n'a ete 
enoncee. Aussi le volume V peut-il etre : 

- un volume fixe dans le referentiel de l'observateur, et done continuellement traverse 
par de nouvelles molecules de fluide. On parle de volume "Eulerien". Dans une 
installation, e'est par exemple le cas d'un recipient, d'une cuve , d'une pompe, 
d'un trongon de tuyauterie. Selon la terminologie de la thermodynamique, e'est 
un systeme ouvert. 

- un volume dont les frontieres bougent au cours du temps. En particulier, si ces 
frontieres bougent avec le fluide, on parle de volume "Lagrangien". Une pro- 
priete fondamentale d'un volume lagrangien est qu'il contient toujours les memes 
particules fluides, meme si celles-ci se deplacent en son interieur. C'est un systeme 
ferme. 

Pour bien differencier les deux, prenons l'image suivante : imaginez une course cy- 
cliste (l'ecoulement), les coureurs etant les particules fluides : 

- le trongon de route entre deux voitures arretees sur le bord est un volume eulerien 
car des cyclistes toujours differents passent continuellement devant ces voitures. 

- le trongon de route (mobile) compris entre la voiture de tete et la voiture-balai 
est un volume lagrangien car il contient toujours les memes cyclistes. II se deplace 
avec les coureurs. 



1.5 Volume fixe ou mobile? 
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On remarquera deux choses importantes : 

- le nombre de cyclistes peut varier dans le volume eulerien (demandez-vous pour- 
quoi) alors qu'il est fixe dans le volume lagrangien. De meme la masse de fluide 
peut varier dans un volume eulerien, mais est constante dans un volume lagran- 
gien. 

- le volume du volume eulerien est constant (les voitures-frontiere sont arretees). 
En revanche lorsque le peloton s'etire, le volume du volume lagrangien augmente. 



La figure 1.3 montre des exemples de volumes euleriens et lagrangiens. 



T 




FlG. 1.3 - En blanc : volume eulerien, en grise : volume lagrangien 



Chapitre 2 



Introduction aux bilans 



2.1 Introduction intuitive 

Le genie des precedes, et en particulier le mouvement des fluides, est regi par des 
principes dits "de conservation" ou encore de "Man" (terminologie chere au genie 
chimique). Un bilan decrit les variations temporelles d'une grandeur au sein d'un 
systeme, que Ton considerera ouvert (au sens de la thermodynamique), c'est-a-dire 
que la grandeur peut rentrer ou sortir du systeme. 

Tout un chacun sait ecrire un bilan, tout au moins lorsque la grandeur fait partie de la 
vie courante. Prenons par exemple le cas ou la grandeur est une quantite d'argent A 
et le systeme mon compte en banque. Pour connaitre la quantite d'argent instantanee 
sur mon compte, j'ecris entre deux instants ti et t 2 : 



Mon compte Mon compte Mes revenus Mes depenses 

en banque - en banque = (ce qui rentre) — (ce qui sort) (2.1) 
en banque a t 2 en banque a t\ entre tiet t 2 entre t±et t 2 

Pour formaliser plus avant ces variations, appelons : 

- G(t) la quantite d'argent instantanee sur mon compte (on supposera a tort que 
c'est une fonction derivable du temps. . .), 

- & e (t) m es revenus instantanes (quantite d'argent qui rentre par unite de temps a 
l'instant t), 

- & s (t) m es depenses instantanees (quantite d'argent qui sort par unite de temps a 
l'instant t). 

Prenons t\ et t 2 infiniment proches soit t\ — t et t 2 — t + dt. Le bilan s'ecrit : 

G(t + dt) - G(t) = <Z> e (t)dt - <Z> s (t)dt 
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2.1 Introduction intuitive 




An 

^ £^ - tjW (2.2) 

Varia^deG 
pendant ac. 



Le membre de gauche est une variation (positive ou negative) qu'on appellera terme 
d'accumulation et les termes de droite sont des termes de flux, respectivement 
entrant et sortant. Les trois grandeurs s'expriment par exemple en Euro / jour. 

Dans cet exemple, il y a une quantite qui se conserve : l'argent (jusqu'a preuve du 
contraire. . .). C'est-a-dire qu'elle ne peut ni disparaitre ni apparaitre, elle ne fait que 
s'accumuler dans le systeme ou entrer / sortir du systeme. En physique c'est le cas 
de l'energie et de la masse. 

Prenons un autre exemple : le systeme est un pays (Fig. |2.1[ ), avec ses frontieres, la 
grandeur est le nombre de personnes instantane G(t) dans ce pays. Ses variations 
peuvent provenir de 3 causes : 

- des personnes entrent ou sortent du pays en traversant ses frontieres, 

- des personnes naissent, 

- des personnes meurent. 

On voit que les deux dernieres causes n'existaient pas dans le cas de l'argent. Ici il 
peut y avoir creation ou destruction de la grandeur etudiee. Le bilan s'ecrit alors, 
en reprenant les memes notations pour les flux entrant et sortant, et en notant 
respectivement R + (t) et R~(t) les nombres de naissances et de deces par unite de 
temps : 



2.2 Transport diffusif et convectif 
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dG 
~dt 



v ■ , ■ An Flux de G Flux de G 

pendZdt entrant - SOTtant - (2-3) 

+ R+{$ - R~{t) 
Production de G Destruction de G 



Les termes R + et R sont classiquement appeles terme source et terme puit. 

En physique, des grandeurs telles que G peuvent etre par exemple un nombre de 
moles (par reaction chimique), la quantite de mouvement (a cause des forces exte- 
rieures), l'entropie (selon le second principe). Ces quantites ne se conservent pas a 
proprement parler puisqu'on peut les creer ou les detruire Q mais on parle malgre 
tout d'equation de conservation. 



Notons que la distinction entre termes sources/puits et termes de flux est claire du point de 
vue physique, elle peut etre arbitraire, mathematiquement parlant. Nous verrons par exemple 
que le premier principe de la thermodynamique, qui ne fait qu'exprimer la conservation de 
l'energie peut etre enonce en prenant pour grandeur energie interne + energie cinetique et en 
prcnant pour termes sources/puits la puissance des forces exterieure + puissance thermique, 
ce qui conduit a une equation du type (2.3). Mais la puissance thermique elle-meme peut etre 
un flux commc nous le verrons ulterieuremcnt. 



2.2 Transport diffusif et convectif 



Nous avons vu qu'un volume peut se vider ou se remplir par ses frontieres d'une 
grandeur extensive G : d'argent (flux financier), de personnes (flux migratoire) . . . 
Dans un fluide, les grandeurs extensives sont la masse, la quantite de mouvement 
ou l'energie (cinetique, interne). Lorsque ces grandeurs traversent la frontiere du 
volume, on parle de transport. 

De fagon generate, il convient de bien distinguer deux types de transport : 

- le transport diffusif intervient meme en l'absence de mouvement visible a 
1'eclielle macroscopique. 

- le transport convectif lie directement au mouvement macroscopique du fluide 



2.2.1 Qu'est ce qu'un flux? 



II convient tout d'abord de mesurer la "vitesse" a laquelle le volume se remplit ou 
bien se vide. Par exemple pour un flux migratoire on pourra parler de nombre de 

: dans le cas de l'entropie, le second principe de la thermodynamique en interdit la destruction. . . 
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2.2 Transport diffusif et convectif 



personnes qui entrent ou qui sortent du pays par annee, ou par toute autre unite 
de temps. Pour un fluide on cherchera par exemple a connaitre la quantite d'eau 
qui coule d'un robinet par unite de temps. Ce type de grandeur s'appelle un flux. 
II est homogene a l'unite de la grandeur G par unite de temps : des kg/s pour un 
flux massique, des Joule/s (ou Watt) pour un flux d'energie etc. C'est un vecteur 
puisque le transport va d'un point vers un autre. 

En general la grandeur G ne rentre (ou sort) pas partout a la meme vitesse du 
volume V. Pour cela il est interessant de ramener le flux a une unite de surface. On 
parle alors de densite de flux (qui est aussi un vecteur). 

2.2.2 Transport diffusif (parenthese) 

Ce type de transport est provoque par le desequilibre d'une grandeur au sein du 
milieu, le systeme tentant de revenir a l'equilibre. Tout le monde sait par exemple 
que lorsqu'il existe un point chaud et un point froid dans un corps, la chaleur est 
transported du point chaud vers le point froid pour tenter d'homogeneiser les tempe- 
ratures : il s'agit d'un flux d'energie, et l'on parle de diffusion thermique. Ce flux 
d'energie n'est pas transports par un mouvement global de fluide, mais aux mouve- 
ments a l'echelle microscopique. La densite de flux energetique correspondante est 
donne par la loi de Fourier (cf. cours de transfert de chaleur) : 

= -AgradT 

Le parametre A est appele conductivity thermique et son unite S.I. est le W/m/K. 

De meme, si l'on verse du lait dans du cafe tres lentement pour eviter tout mou- 
vement, on constate que les deux finissent par se melanger parfaitement : dans un 
melange de deux especes A et B, l'espece A diffuse des regions ou elle est la plus 
concentree vers celle ou elle Test le moins. Ici encore, ce flux de matiere n'est pas 
transports par un mouvement global de fluide, mais est lie uniquement aux collisions 
entre les differentes molecules a l'echelle microscopique. La densite de flux massique 
de l'espece A est fournie par la loi de Fick (cf. cours de transfert de matiere). 

Ja = -D A bP grades 

ou uja est la fraction massique de A, et Dab est appele coefficient de diffusion de A 
dans B (en m 2 /s). 

La quantite de mouvement dans un fluide en mouvement diffuse egalement, a cause 
du frottement des couches fluides les unes sur les autres, phenomene appele visco- 
site. Nous y reviendrons lors de l'etude des forces exercees sur un fluide. 

Les transferts diffusifs, et notamment les deux formules ci-dessus, ne nous serviront 



2.2 Transport diffusif et convectif 
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pas dans le cadre de ce cours car : 

- les transferts de matieres diffusifs n'interviendront pas parce qu'ils ne concernent 
que des melanges de plusieurs fluides, non considered dans ce cours. 

- les transferts de chaleur diffusifs interviendront comme des termes source/puit de 
chaleur (Q) dans l'expression du premier principe, mais ne seront pas explicites. 

- les transferts diffusifs de quantite de mouvement interviendront implicitement 
dans les forces appliquees au fluide, considerees comme des sources/puits de quan- 
tite de mouvement. 



2.2.3 Flux convectif 



Le flux convectif (du latin cum = avec, vectare = transporter) est lie au transport 
d'une quantite par le mouvement du fluide. Tout le monde a un jour ou l'autre 
remue l'eau dans une baignoire pour transporter l'eau chaude coulant du robinet 
vers l'autre bout de la baignoire. L'agent de transport est dans ce cas le fluide lui- 
meme et est bien plus rapide que le flux diffusif. C'est pour la meme raison que 
Ton remue son cafe pour melanger le sucre ou le lait, sa soupe ou une sauce dans 
une casserole pour la refroidir, que Ton agite les reacteurs chimiques. Rappelez-vous 
egalement combien il est difficile de melanger deux peintures de couleurs differentes : 
s'il fallait attendre qu'elle se melangent par diffusion, il faudrait attendre des annees, 
c'est pourquoi il faut remuer le melange. 

Le but de cette section est de quantifier ce flux convectif pour n'importe quelle 
grandeur traversant la surface frontiere S d'un volume V eulerien, c'est-a-dire 
fixe dans l'espace. 

Prenons pour commencer l'exemple simple suivant : de la fumee sort d'une piece 
fermee par une fenetre a cause d'un vent de vitesse v(t), uniforme et perpendiculaire 



a la fenetre (Fig. 2.2), par exemple cause par un ventilateur. Intuitivement, on 
congoit bien qu'il sortira, pendant un temps donne, d'autant plus de fumee que 
celle-ci est concentree d'une part, et que le vent est violent d'autre part. 

Passons a revaluation de la masse de fumee dG sortant pendant un temps dt : c'est 
celle qui se trouve dans un volume S dl ou dl = v dt est le trajet parcouru par la 
fumee pendant un temps dt (Fig. 2.2[ ), soit, si nous appelons g la densite de la fumee 



(la masse de fumee par unite de volume) 

= 9 x x x ^jttji 

k S kg/m 3 m2 m / s s 

Imaginons maintenant que la densite de la fumee g et/ou la vitesse v ne soit plus 



uniforme sur toute la surface de la fenetre (Fig. 2.3), mais depende de la position 
(x,y) du point de la fenetre considere : il faut alors calculer les quantites infinite- 
simales rentrant par chaque petit element de surface dS = dx dy et integrer sur 
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2.2 Transport diffusif et convectif 




v(t) 




dl=vdt 



Fig. 2.2 - Fumee sortant d'une piece par une fenetre. 

toute la surface. La quantite de fumee d 2 G sortant de la piece q par dS est la fumee 
contenue dans le volume dS x v dt, soit : 



d 2 G = gdS x vdt 




v{x,y,t) 




v(x, y)dt 



Fig. 2.3 - Fumee sortant d'une piece par une fenetre avec une vitesse non homogene. 
soit, en integrant sur tout la surface de la fenetre : 

dG = dt x J J gv dS 
' s 



2 On le note ainsi car dS et dt etant tous deux infiniment petits, d 2 G est un infiniment petit 
d'ordre 2. 



2.2 Transport diffusif et convectif 
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X 



Fig. 2.4 - Fumee sortant d'une piece par une fenetre avec une vitesse non homogene 
et non perpendiculaire a la fenetre. 



Notons que cette expression, sous reserve que v represente la composante (positive 
ou negative) du vecteur v suivant l'axe sortant de la piece, autorise qu'en certains 
points de la fenetre, la fumee re-rentre dans la piece : l'integrale ci-dessus est alors 
le bilan entre la masse sortante et la masse entrante dans la piece. 

II reste a voir le cas le plus general ou la vitesse n'est pas orthogonale a la surface 



(Fig. 2.4 ). Considerons un petit element dS de la fenetre et designons par n le vecteur 
normal a dS, oriente de I'interieur vers Vexterieur de la piece (on parlera de vecteur 
normal sortant). 

par dS est la quantite de fumee dans le volume pointille sur la figure. Ce volume 
vaut dS x (vdi.n) et done : 



d 2 G = gdS x vrft.n 



et le flux sortant de la piece pendant dt est obtenu en integrant sur tous les elements 
de surface dS : 



dG = dtx J J gv.n dS (2.4) 

s 

Rappelons que les grandeurs a I'interieur de l'integrale dependent du point de la 
surface consideree : g parce que la concentration de fumee peut varier d'un point a 
un autre, v parce que la vitesse du fluide peut varier d'un point a un autre, et n 
parce que la surface peut etre courbe, auquel cas le vecteur normal varie d'un point 
a un autre. 
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2.2 Transport diffusif et convectif 



Notons ensuite que cette definition fonctionnerait aussi si la fumee entrait dans la 
piece au lieu d'en sortir : on aurait alors v.n < en tout point, et la quantite de 
fumee sortant de la piece dG est aussi negative. Dans le cas general oil la fumee peut 
entrer en certains points (v.n < 0) et sortir en d'autres (v.n > 0), on voit que la 



formule (2.4) fournit le bilan entre la fumee sortante et la fumee entrante 



pendant l'intervalle de temps dt. C'est done tres precisement l'oppose de la 



grandeur dt x [$ e (i) — $ s (^)] intervenant dans l'equation de bilan (2.3), et nous 
avons done : 



Par unite de temps, un fluide de champ de vitesse v transporte a l'interieur d'un 
volume V fixe dans l'espace une quantite nette de G egale a : 



$ e (0 - = - J J gv.ndS (2.5) 

s 

oil g est la quantite de G par unite de volume. 



Dans le cas oil le volume V verrait ses frontieres bouger, au lieu d'etre fixes, il est 
facile de generaliser ce resultat en remplagant le champ de vitesse v par v — w, 
ou w est la vitesse locale en chaque point de la frontiere. Dans le cas extreme ou 
w = v, e'est-a-dire quand la frontiere bouge avec le fluide (V est alors un volume 
Lagrangien), on voit clairement que le flux convectif est nul : le fluide ne peut 
transporter une grandeur a travers la frontiere puisque vu depuis celle-ci, il est 
immobile. 



2.3 Bilan d'une grandeur volumique dans un milieu continu 
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2.3 Bilan d'une grandeur volumique dans un mi- 
lieu continu 

2.3.1 Ecriture generale 



II suffit maintenant d'introduire ce terme de flux dans le bilan (2.3). 



Les variations d'une grandeur volumique g dans un volume V limite par une 



surface S (Fig. 2.5) s'ecrit : 



dG 
~dl 



d 
dl 




gdV 



v 



Accumulation de 
G 



gv.n dS 



s 



Flux convectif de 
G 

entrant - sortant 



5-3 

Production de G 



Destruction de G 



(2.6) 




Fig. 2.5 - Bilan d'une grandeur g sur un volume V 

Rappelons que cette formule est liee a l'orientation du vecteur normal a la frontiere 
vers l'exterieur. Si le vecteur normal etait pris vers l'interieur nous aurons un signe 
+ devant l'integrale de surface. 



Les trois equations de bilan regissant le mouvement des fluides sont toutes de ce type, 
respectivement avec g = p pour la masse, g = pv pour la quantite de mouvement 
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2.3 Bilan d'une grandeur volumique dans un milieu continu 



et g = p(e + v 2 /2) pour l'energie. Pour chacune d'entre elles, il nous reste juste a 
specifier quels sont les termes source et puit. Ces derniers resultent de trois grands 
principes physiques : 

- l'impossibilite de creer ou detruire de la matiere (conservation de la masse), 

- la loi de la dynamique de Newton (conservation de la quantite de mouvement), 

- le premier principe de la thermodynami que (conservation de l'energie). 



2.3.2 Une geometrie particuliere : le tube de courant 



Un tube de courant (Fig. 2.6) est un tube fictif dont les parois laterales S\ ai sont 
en tout point tangentes au vecteur vitesse et comportant une section d'entree S e et 
une section de sortie S s , toutes deux droites. II peut s'agir par exemple d'un tube 
reel car la condition d'etancheite impose que la composante normale de la vitesse 



sur une paroi solide (immobile) soit nulle (cf. section 4.2) 



n 




Fig. 2.6 - Tube de courant. 



Dans ces conditions l'integrale de surface se limite aux sections d'entree et de sortie, 
car v.n est nul partout sur la surface laterale Si at , et le bilan d'une grandeur g s'ecrit, 
en decomposant la surface S en S e + S s + S\ at '■ 



dG 

Accumulation de 
G 



gv.n dS — 



Se 



G e = Quantite 

> de G 
rentrant dans le 
tube par S e par s 



gv.n dS 



Ss 



G s = Quantite 

> de G 
rentrant dans le 
tube par S s par s 



Taux de 
production de G 



R-(t) 



Taux de 
destruction de G 



(2.7) 



2.3 Bilan d'une grandeur volumique dans un milieu continu 
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La premiere integrate de surface est negative car v.n < partout sur S e et la seconde 
est positive car v.n > partout sur S s . On peut done les exprimer en fonction des 
deux grandeurs G e et G s , toutes deux positives representant des debits de G (masse, 
energie,. . .) respectivement entrants et sortant par le tube. 

Le tube de courant est une notion tres importante car il permet de representer, au 
sein d'un reseau de fluide, n'importe quel systeme ayant une entree et une sortie 
(troncon de tube, pompe, vanne, coude . . .). 

Un cas particulier tres frequent est le regime permanent ou les proprietes de 
l'ecoulement ne dependent plus du temps. C'est le cas de nombreuses installations 
industrielles, du circuit de freon de votre refrigerateur, de votre cuisiniere a gaz 
etc. . .. Dans le cas du tube de courant on obtient une expression particulierement 
simple du bilan : 



= G^ + R + (t) - R~(t) 

Le G qui sort Le G qui rentre Le G qui est Le G qui est 

produit detruit 

(2.8) 



2.3.3 Une approximation utile : l'ecoulement piston 



Lorsque Ton traite des ecoulement dans des tubes reels, qui sont un cas particulier 
du tube de courant, il est tentant de supposer que le champ de vitesse ne varie 
pas transversalement dans les parties du tube ou la section est constante : on parle 
alors d'ecoulement piston car tout se passe comme si chaque tranche de fluide se 



deplagait comme un piston solide (Fig. 2.7). 



Se 




1 ► 


► v 


n 


► 

► 



5., 



n 



Fig. 2.7 - Modele de l'ecoulement piston. Le vecteur vitesse est uniforme dans toute 
section du tube. 



Dans ce cas, les integrates sur les surfaces d'entree et de sortie intervenant dans le 



bilan (2.7) se simplifient car le produit scalaire v.n est maintenant constant sur S e 



ou S s et peut etre sorti de l'integrale. Le debit de sortie s'ecrit par exemple : 
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2.3 Bilan d'une grandeur volumique dans un milieu continu 



Gs = JJ ^ n dS = Vs J J 3 dS = gsVsSs ( 2 ' 9 ) 

Ss Ss 

ou v s est le module de la vitesse a la sortie du tube, et g s est la valeur moyenne de 
g sur la section S s . On a de meme : 

a. = -// 9 v.ndff = «.//,dff = fc « A (2.10) 

S e S e 



Si un ecoulement dans un tube peut etre considere comme un ecoulement piston, 
nous pouvons ecrire, en tout point du tube le debit d'une grandeur g sous la 
forme : 

G = gvS (2.11) 

oil v est la vitesse du fluide au point considere, S la section du tube et g la 
valeur local de la grandeur volumique considered. 



Dans ces conditions, le bilan (2.7) entre l'entree et la sortie du systeme prend une 
forme particulierement simple, puisque les integrates ont disparu, au profit de values 
moyennees sur les sections. On obtient 



Le bilan d'une grandeur g sur un tube de courant dans l'hypothese d'un ecou- 
lement piston s'ecrit : 

dG 

— = geVeSe-gsVsSs+R + (t) - R~ {t) (2.12) 

G e G s 

Dans cette equation v e et v s sont des vitesses moyennes sur les sections d'entree 
et de sortie. 



II reste a savoir si une telle approximation est justifiee. En fait, un profil de vitesse 
uniforme dans la section du tube est impossible : en effet, comme nous le verrons 
plus tard, le frottement visqueux du fluide sur la paroi impose que le champ de 
vitesses est necessairement nul sur celle-ci. Le champ de vitesse doit done en realite 
passer continument de a la paroi a sa valeur maximale au centre du tube. Selon 
que l'ecoulement est laminaire ou turbulent, le profil de vitesse est soit parabolique 



(on parle alors d'ecoulement de Poiseuillepl) , soit aplati (Fig. 2.8). 



3 Ce resultat est rclativement simple a demontrer a partir des equations de Navier-Stokes, qui 
releve de la seconde partie du cours consacree aux fluides reels (c.a.d pour lesquels on prend en 
compte le frottement visqueux). 



2.3 Bilan d'une grandeur volumique dans un milieu continu 
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Fig. 2.8 - Profils de vitesse dans un tube pour un ecoulement laminaire (a gauche), 
et turbulent (a droite). 

Cependant, il est toujours possible de representer ces types d'ecoulement reels par 
des ecoulements piston dont la vitesse est une moyenne de la vitesse sur la section 
du tube. On pose alors : 



III 



, , v.j i dS 

section 



et Ton considere que la vitesse prend cette valeur en tout point du tube. Pour un 
ecoulement laminaire, la theorie de Poiseuille montre que cette vitesse moyenne est 
egale a 1/2 de la vitesse max, realisee au centre du tube. 



2.4 Equations de conservation pour un fluide 
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2.4 Equations de conservation pour un fluide 
2.4.1 Conservation de la masse 

Le principe physique stipule que (hormis pour des reactions nucleaires), la masse de 
matiere se conserve. Par consequent un volume V fixe dans l'espace peut accumu- 
ler de la matiere ou en echanger avec l'exterieur mais pas en creer ni en detruire. 
L'equation de bilan de masse d'un fluide ne comporte ni terme source, ni terme puit. 

Utilisant les resultats du chapitre precedent, nous avons done, en prenant G = M 
la masse du volume V et g = p la masse volumique : 




2.4.2 Conservation de la quantite de mouvement 
2.4.2.1 Rappel de mecanique 

On peut questionner la pertinence du terme "conservation". II convient cependant 
de se rememorer comment est definie la notion de force en mecanique : une force 
est associee a l'echange de quantite de mouvement entre deux systemes. Si ces deux 
systemes sont isoles de l'exterieur, la somme de leurs quantites de mouvement reste 
constante au cours du temps, quelle que soit leur interaction mutuelle. Par exemple 
un couple de patineurs sur la glace est quasiment isole de l'exterieur si Ton peut 
negliger le frottement des patins. S'ils se tiennent par la main, ils interagissent et 
echangent de la quantite de mouvement, mais la somme des deux reste constante. De 
meme deux boules de billard qui se choquent echangent de la quantite de mouvement 
lors du choc, mais la somme des deux reste constante. 

C'est en ce sens qu'il faut comprendre l'expression "conservation de la quantite de 
mouvement" : on peut dire que le fluide subit (ou exerce) une force exterieure, mais 
aussi qu'il echange de la quantite de mouvement avec l'exterieur, de telle sorte que 
la somme des quantites de mouvement fluide + exterieur reste constante. II est 
important de retenir cette notion pour bien comprendre la physique sous-jacente 
aux equations qui vont suivre. 

On retiendra done : 

Une force subie par le fluide correspond a un accroissement de sa quantite de 
mouvement. Une force exercee par le fluide correspond a une diminution de sa 
quantite de mouvement 
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2.4 Equations de conservation pour un fluide 



2.4.2.2 L'equation generale 



On considere toujours le volume V fixe dans l'espace. Les termes sources R + — R 
sont les forces exterieures ^F cxt , et l'equation de bilan (2.6) s'ecrit, avec G — P, 
et g = pv : 




Cette relation appelle plusieurs remarques : 

1. Le lecteur ayant des notions en mecanique s'etonnera de la difference entre 
cette equation et la loi de la dynamique bien connue dP /dt — Yl F ex t : il y 
a un terme de flux en plus ! Ce terme de flux vient du fait que contrairement 
aux systemes etudies habituellement, V est un systeme ouvert (ou eule- 
rien) qui ne contient jamais les memes particules materielles ! Aussi 
les particules fluides rentrant et sortant du volume V transportent-elles avec 
elles leur quantite de mouvement. C'est l'origine du terme de flux. Si nous 
avions pris un volume V(t) lagrangien, qui suit le fluide dans son mouvement, 
et contient done toujours la "meme matiere", les termes de flux disparaitraient 
(pourquoi ?) et on retrouverait bien la loi de la dynamique dans sa formulation 
classique. 

2. C'est une equation vectorielle et ce sont done en fait trois equations faisant 
intervenir les trois composantes de la vitesse. 

3. Le terme pv(v.n) doit etre lu comme le "produit scalaire de v par n, ensuite 
multiplie par le vecteur pv. 



2.4.3 Conservation de l'energie 

2.4.3.1 Rappel : premier principe en systeme ferme. 



Le premier principe de la thermodynamique stipule que, pour un systeme ferme : 



A(U + K) = W + Q (2.15) 

oil U energie interne, K energie cinetique, W travail des forces exterieures et Q 
la quantite de chaleur echangee (comptee positivement si le systeme regoit de la 
chaleur). Toutes ces grandeurs se comptent en Joule. Ecrivons le premier principe 
pendant un intervalle de temps dt, nous obtenons 
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-(U + K) = W + Q (2.16) 

ou W et Q sont respectivement la puissance des forces exterieures et la puissance 
thermique echangee, comptees en Watt (attention la notation avec un point au- 
dessus ne signifie pas qu'il s'agisse d'une derivee temporelle). 



2.4.3.2 Adaptation au systeme ouvert. 



Nous ne pouvons appliquer directement ce principe a notre volume V puisqu'il est 
traverse continuellement par de la matiere differente : c'est un systeme ouvert. Cette 
matiere entrante et sortante transporte de l'energie cinetique et interne vers ou hors 
du volume V. Comment appliquer le premier principe a u n syst eme ouvert ? Tout 
simplement en remarquant que l'enonce en systeme ferme ( 2.16[ ) indique que W et 
Q sont des termes sources/puits pour les variations de la grandeur "energie interne 
+ energie cinetique". Pour exprimer le premier principe en systeme ouvert, il sufEt 
de rajouter les termes de flux convectif dans l'equation de bilan generique (2.6) avec 
G = U + K, g = p(u + v 2 /2) et R+ - R~ = W + Q, soit : 




Ici comme pour la quantite de mouvement, on remarquera que si le volume V est 
un volume materiel, qui suit le fluide dans son mouvement (et est done un systeme 
ferme) le terme de flux disparait, et l'on retrouve l'expression classique du premier 
principe en systeme ferme. 



2.4.4 Synthese 



Rappelons les equations de bilan obtenues pour la masse, la quantite de mouvement 
et l'energie, ecrites toutes les trois sous la forme d'un terme d'accumulation, un 
terme de flux sortant (algebrique, done negatif pour un flux entrant) et un terme 
source (egalement algebrique, c.a.d positif si production, negatif si destruction) : 
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2.4 Equations de conservation pour un fluide 



Accumulation + Flux sortant = Source 

-Flux entrant (2.19) 



d 
dt 




d_ 
dt 




pdV + J J pv.n dS = (2.20) 

r s 

pvdV + [[(pv)v.ndS = ^F cxt (2.21) 




V s 



v 2 \ , f f ( v 2 



( p lu+-\ dV + pl u +-)v.ndS = W + Q (2.22) 



v 



Pour detailler plus avant les equations de conservation dans un fluide, il nous manque 
l'expression des forces exterieures F cxt agissant sur le fluide, ainsi que la puissance 
thermique Q. 



Le prochain chapitre decrit l'ensemble des forces (les plus classiques) exercees sur le 
fluide. 



Chapitre 3 



Forces exercees sur un fluide 



3.1 Introduction 



Nous considerons toujours un volume de fluide V eulerien. 

II est d'usage de distinguer deux types de forces dans les milieux continus : 

- les forces volumiques qui s'exercent sur chaque particule fluide interne au vo- 
lume V. 

- les forces surfaciques ou de contact qui s'exercent sur la frontiere S du vo- 



3.2 Force volumiques 

Pour un fluide, en voici une liste possible : 

- le poids 

- les pseudo-forces d'inertie d'entrainement et de Coriolis en referentiel non-galileen 
(importantes pour les mouvements de l'atmosphere par exemple). 

- les forces electriques et magnetiques pour les fluides charges (suspensions de par- 
ticules magnetiques, magma) [j] 

Nous nous limiterons ici au poids. Exprimons le poids du volume de fluide V : chaque 
volume elementaire dV a pour poids pgdV de telle sorte que le poids du volume V 
s'exprime par : 



lume V. 





v 



^^Cette disciplique complexe porte le nom de "magnetohydrodynamique' 
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3.3 Forces de contact : pression 



On notera que p et g ne peuvent etre sortis de l'integrale dans le cas le plus general : 
p parce que le fluide peut presenter des variations de masse volumique, et g parce 
qu'a l'echelle d'un ocean ou de l'atmosphere, g ne garde pas une direction constante. 



3.3 Forces de contact : pression 
3.3.1 Origine microscopique. 

Nous avons mentionne en introduction le caractere granulaire des fluides a l'echelle 
microscopique, et le fait que ces "grains" sont en perpetuel mouvement les uns par 
rapport aux autres. 

Supposons que le fluide est limite par une paroi solide (par exemple celle d'un re- 
cipient). A chaque instant, certaines des molecules du fluide, dans leur mouvement 
desordonne, vont aller taper et rebondir sur cette paroi (imaginez une multitude de 
balles lancees aleatoirement dans un gymnase) (Fig. |3.1|). 



\ \ \ I / 




Fig. 3.1 - Choc de molecules d'un fluide sur une paroi solide. 
Lors de ce rebond, conformement aux lois de Newton, la molecule echange de la 
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quantite de mouvement avec la paroi, et exerce done sur cette derniere une force 
proportionnelle a {sa quantite de mouvement avant le choc /xvi} - {celle apres le 



choc /iv'jj. On voit sur la figure 3.2 que cette force est perpendiculaire a la paroi et 
dirigee vers son interieur. 




F OC /LtVi 



Fig. 3.2 - Echange de quantite de mouvement lors du choc d'une molecule de fluide 
de masse /i sur une paroi solide. 

La force moyenne qu'exerce le fluide sur la paroi, obtenue en comptant tous les chocs 
de molecules par unite de temp^] est done egalement normale a la paroi et dirigee 
vers son interieur. Cette force est une force de pression. 

II est bien evident que le nombre de chocs subi est proportionnel a la surface de la 
paroi. Si Ton considere une surface deux fois plus petite, le nombre de chocs, et done 
la force resultante, seront egalement deux fois plus petits. La force de pression est 
done proportionnelle a la surface considered. II est done commode dans ces conditions 
de definir une force par unite de surface appelee pression (l'unite S.I. est le Pascal : 
1 Pa = 1 N/m 2 = 1 kg.m _1 .s -2 ). Comme l'orientation de cette force est parfaitement 
definie (normale et rentrante dans le systeme subissant la force de la part du fluide) , 
un scalaire p suflit pour la definir. 

Quotidiennement, notre corps subit les memes chocs que la malheureuse paroi de la 



figure 3.1 de la part des molecules d'air environnant. La pression correspondante, 
dite pression atmospherique est egale a 101325 Pa, ce qui correspond grossomodo a 
une masse de 1 kg/ cm 2 ! 



Sur les figures |3.1| et |3.2[ rien n'empeche de considerer autre chose qu'une paroi 
solide : on peut aussi bien considerer un volume de fluide identique ou non a celui 
qui exerce la force. Les molecules du fluide de gauche exerceront les memes chocs sur 
le fluide de droite, et done la meme pression. II est done possible de la pression en 
un point dans un fluide. Nous allons eclaircir cette notion dans la section suivante. 



2 Un calcul approche simple peut etre effectue pour un gaz parfait, et fournit l'equation d'etat 
des gaz parfaits. 
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3.3 Forces de contact : pression 



3.3.2 Equilibre d'une colonne d'eau. 



Imaginons 3 cylindres de liquide immobiles empiles les uns sur les autres (Fig. 3.3), 



que nous nommerons 1, 2 et 3. Le dessus du cylindre 1 est soumis a la pression 
atmospherique p a tm- Nous allons considerer l'equilibre du cylindre du milieu, pour 
l'instant seulement dans la direction verticale. II est soumis a son propre poids, a une 
force de pression Fi/ 2 de la part du cylindre 1 et d'une force de pression verticale 
F3/2 de la part du cylindre 3. L'equilibre s'ecrit done : 



0- 



zi = Li - 



z\ = Li + L 2 - 




Fig. 3.3 - Equilibre d'un cylindre de fluide. 



M 2 g + Fi /2 + F3/2 = 



Considerons maintenant l'equilibre du cylindre 1 : il est soumis a son poids, a la 
force de pression F exercee par la pression atmospherique et a la force —Fx/2 de la 
part du cylindre 2, d'apres la loi de Faction et de la reaction. Son equilibre s'ecrit 
done : 



M lg - Fi /2 + F = 



3.3 Forces de contact : pression 
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De ces deux equations, et en remarquant que F 3 / 2 = — F 2 /3, nous deduisons les 
forces F1/2 et F 3/2 : 



F V 2 = F + M lg 

F2/3 = F + (Mi + M 2 )g 

ce qui montre que chaque cylindre subit de la part de son superieur une force egale 
au poids de toute la colonne le surplombant + la force exercee par l'exterieur sur le 
cylindre 1 qui se transmet dans toute la colonne. 

De plus, appelons S la section des cylindres, p la masse volumique du fluide, et 
Li, L 2 , L 3 leurs longueurs respectives. Les deux equations d'equilibres precedentes 
s'ecrivent alors : 



F1/2 

S 
F2/3 



Fo 

S 
Fo 

S 



+ p(L 1 + L 2 ) 



(3.2) 



Mais la force Fo exercee par l'atmosphere est liee aux chocs des molecules d'air sur la 
surface superieure du cylindre 1, et c'est done une force de pression, proportionnelle a 
la surface S et orientee vers le bas. On peut done l'ecrire sous la forme F = Pa,tmSe z . 
Les equations precedentes deviennent : 



F1/2 T 

J (3-3) 

i? 2/3 

-^r~ = Pat m e z + p(Li + L 2 )g 

Ces expressions montrent clairement que si on double ou decuple la section des 
cylindres, les forces changent certes dans le meme rapport, mais pas le rapport F/S !. 
Les forces Fi/ 2 et F 2 /3 sont aussi des forces de pression : elle peuvent ainsi s'exprimer 
comme le produit de la pression p locale par la section de contact. On voit de 



plus d'apres les expressions (3.3) que l'amplitude de la force exercee sur la surface 
superieure d'un cylindre quelconque ne depend que de la profondeur z de cette 
surface (z = L\ pour le cylindre 2, z = L\ + L 2 pour le cylindre 3). La pression 
ne depend done que de la profondeur et Ton peut ecrire les formules ( |3.3[ ) dans une 
formule unique : 



Pi*) = Patm + PgZ 



(3.4) 
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3.3 Forces de contact : pression 



en prenant l'axe z oriente dans le sens de la pesanteur et le zero au niveau de 
la surface. Ce resultat est classique, il indique que la pression augmente avec la 
profondeur dans un fluide immobile. Nous retrouverons ce resultat plus generalement 
ci-dessous. 



3.3.3 Generalisation. 



Nous allons recapituler dans une formule ce que nous avons appris dans les sections 
precedentes : la force de pression exercee sur un systeme est : 

- proportionnelle a la surface 

- normale a la surface 

- orientee vers l'interieur du systeme 

Considerons done un volume V quelconque de fluide immobile entoure du meme 



fluide (Fig. 3.4). Sur chaque element dS de la frontiere de V, le fluide exterieur 
applique une force normale rentrant dans le volume V, definie par la pression p 
locale (dependant done de x, y et z) : 



On definit la force de pression sur un element de surface dS par 

dF = -pn dS (3.5) 
et done, pour la surface S toute entiere : 

F = ff -pn dS (3.6) 



s 



Rappelons que l'unite S.I. pour la pression est le Pascal egal a un N/m 2 . La pression 
atmospherique au niveau du sol vaut 101325 Pa, unite appelee aussi atmosphere. Le 
Pascal est done une unite tres petite. Le bar, unite de pression dans le systeme CGS 
vaut 10 5 Pa (d'ou les fameux 1013 mbar de la meteo), soit environ une atmosphere. 

Le psi (Pound per Square Inch) est egal a une livre par pouce-carre, soit 0.4536 kg 
/ (2.54 x 10~ 2 m) 2 ~ 6897 Pa. On pourra retenir que 1 bar ~ 14.5 psi (verifiez 
lorsque vous gonflerez les pneus de votre voiture :-). On ajoute generalement une 
lettre au symbole psi selon ce que Ton mesure : 

- psia represente une pression absolue (p) 

- psig represente une pression relative par rapport a la pression atmospherique 

(P - Patm) 

- psid represente une pression differentielle entre deux points (p2 — Pi) 



3.3 Forces de contact : pression 
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On trouve egalement des unites liees a des hauteurs equivalentes de liquide (voir 
ci-dessous). Le torr, ou mm de mercure est la pression exercee par une colonne d'un 
millimetre de mercure. II vaut 133.32 Pa, et on retiendra que 1 atmosphere vaut 760 
mm Hg ou 760 torr. Le torr est bien adapte aux systemes sous vide. 



3.3.4 Loi de l'hydrostatique. 



Nous connaissons maintenant les deux forces exercees sur un fluide immobile : le 
poids et les forces de pression. Done : 



L'equilibre mecanique d'un volume V de fluide immobile s'ecrit : 

pgdV + J J -pndS = (3.7) 

r S 




La figure [3T4 illustre cette equation : pour que la resultante des forces de pression 
equilibre le poids, la pression doit etre plus elevee dans les parties basses que dans les 
parties hautes, ce que nous avions demontre pour la colonne de liquide : la pression 
augmente avec la profondeur. 




Fig. 3.4 - Equilibre d'un volume de fluide : la resultante des forces de pression 
contrebalance le poids. 



Pour mieux quantifier ce resultat, utilisons la formule du gradient (B.2) pour trans- 
former l'integrale de surface en integrale de volume. Nous obtenons : 
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3.3 Forces de contact : pression 




pgdV + JJJ(- gradp) dV = 

v v 

qui doit etre vraie pour tout volume V en equilibre, et par consequent 



En tout point d'un fluide immobile, on a : 




pg — gradp = 


(3.8) 


appelee loi de l'hydrostatique. 





line propriete du gradient bien connue est qu'il est orthogonal aux surfaces equiva- 
leurs, c'est-a-dire que f(x,y, z) est constant lorsque Ton se deplace perpendiculaire- 
ment a grad /. Nous voyons done que : 



Les surfaces isobares dans un fluide au repos sont perpendiculaires a g. 



Ce sont done des surfaces spheriques a l'echelle de la terre (si on oublie les defauts 
locaux de g) et des plans horizontaux a l'echelle humaine. Nous allons specifier deux 
applications pratiques. 



3.3.5 Applications. 

3.3.5.1 Variation de la pression avec la profondeur de liquide 

Prenons un axe Oz, oriente vers le bas (dans le sens de g), et l'origine a la surface 



du liquide. L'equation (3.8) se decompose alors en trois equations : 

ox 

OX 

dp 

dz =P9 

et nous obtenons alors, en supposant le liquide incompressible (p = C te ) 



P(z) = Patm + pgZ 
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On retrouve le resultat obtenu pour une colonne de fluide. Avec p atm = 101300 Pa, 
et p = 998 kg/m 3 , on voit que la pression augmente environ de p atm tous les 10 
metres. 



3.3.5.2 Pression dans l'atmosphere 

Prenons maintenant l'axe Oz vers le haut, et l'origine au niveau de la mer. L'equation 
de l'hydrostatique s'ecrit maintenant : 



dz = ~ p{z)g 

ou nous avons ecrit p(z) car Pair est compressible et les variations de p ne peuvent 
etre negligees dans l'integration. Celles sont donnees par la loi des gaz parfaits : 



R m 

ou M a = 28.8 g/mol est la masse molaire de l'air, et R = 8.32 J.mol _1 .K -1 . Si 
nous considerons T = T constant avec l'altitude, on peut integrer l'equation de 
l'hydrostatique en : 



M a g 



En realite, jusqu'a 11000 m d'altitude (troposphere), la temperature diminue avec l'altitude, 
selon une loi approximativcmcnt lincaire (si vous avez deja pris l'avion, vous avez probablement 
ete surpris en entendant le pilote annoncer une temperature exterieure de l'ordre de -40 °C. . .). 
Ainsi, nous avons : 

T(z) = T Q -Bz 

avec B ~ 6.5 x l0~ 3 K/m. L'utilisation de cette loi fournit alors, apres integration : 

/ Bz \M a9 /B.B 



3.3.6 Extension en referentiel non galileen. 

Examinons le cas ou le fluide est immobile dans un referentiel non galileen que 
nous nommerons TV . Cela se produit notamment dans des cuves tournantes (miroir 
liquide), ou par exemple dans un vehicule en acceleration. 

Puisque le liquide est immobile dans 71', il ne subit que la force d'inertie d'entraine- 
ment, celle de Coriolis etant nulle (cf. annexe O). La force d'inertie d'entrainement 
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sur un element de volume dV s'ecrit — p dVa. e et il est facile de voir que la loi de 
l'hydrostatique se generalise en : 



p(g - a e ) - gradp = (3.9) 

Nous voyons done que : 



Dans un referentiel non galileen TZ', d'acceleration a e par rapport a un referen- 
tiel galileen, la loi de l'hydrostatique s'ecrit de fagon identique, a condition de 
remplacer g par une pesanteur "apparente" g' = g — a e . 



Les surfaces isobares deviennent normales a g'. Ce seront par exemple des plans 
inclines dans une voiture en acceleration constante, et des paraboles dans les re- 
ferentiels tournants. Ce resultat est mis a profit pour construire naturellement des 
miroirs liquides parfaitement paraboliques a partir de cuves de mercure en rotation. 



3.3.7 Poussee d'Archimede. 



Nous avons evalue la force de pression resultante sur un volume V de fluide immobile, 
montre qu'elle etait dirigee vers le haut et que puisqu'il y a equilibre, elle est 
l'opposee du poids du fluide. 

Que se passe-t-il si Ton remplace ce volume fluide V par un corps etranger (bulle, 
goutte, corps solide) ? Tant que le fluide est immobile, le champ de pression dans le 
fluide reste le meme (p = p atm + pgz). Les forces de pression sur la surface S sont 



done les memes (celles de la figure 3.4), que V contienne du fluide ou toute autre 
matiere et leur resultante est done l'opposee du poids du fluide qui remplirait ce 
volume V ! D'ou le celebre principe : 
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Tout corps plonge dans un liquide subit une poussee de bas en haut egale au 
poids du liquide deplace (c'est-a-dire du fluide qui remplacerait le corps si ce 
dernier etait absent). 

Dans le cas le plus courant ou la masse volumique du fluide p peut etre considered 
constante (et si g uniforme), on obtient : 

F a = -pVg (3.10) 

Cette force est appelee force ou poussee d'Archimede^ On retiendra qu'elle 
est toujours opposee au champ de pesanteur g (apparent g' dans le cas d'un 
referentiel non galileen). 



"Archimcde (287 av. JC - 212 av. JC) : sans doute le plus grand scientifique de l'antiquite. 
Outre le fameux principe, qu'il aurait decouvert en prenant son bain (et le Eureka consecu- 
tif. . .), on lui doit de nombreuses inventions comme le lcvier, ou la vis sans fin (qui porte 
son nom). II inventa de nombreuses machines de guerre pour defendre Syracuse. II fut egale- 
ment prolifique dans le domaine des mathematiques : on lui doit notamment la spirale et des 
estimations assez precises de 7r. II fut tue par un soldat romain qui ignorait qui il etait. 



Applications : le corps etranger occupant V est soumis egalement a son poids p s Vg, 
ou p s est la masse volumique du corps. La resultante de son poids et de la poussee 
d'Archimede est done : 



F a + V = (p s -p)Vg 
= {d-l) P Vg 

ou d = p s / p est un nombre sans dimension, qu'on appelle densite lorsque le fluide 
est de l'eau (p = 998 kg.m~ 3 ). La force nette est done vers le haut si le corps est 
plus leger que le fluide (d < 1), et vers le bas s'il est plus lourd (d > 1). On peut 
aussi interpreter cette formule comme si le corps avait une masse modifiee (d — l)pV. 
C'est ainsi que Ton se sent "plus leger" dans l'eau (le corps humain a une densite 
moyenne de l'ordre de 1.07) 

La plupart des bois sont plus legers que l'eau et done y flottent. L'ebene (d=l.l) et le chene 
rouge d'Australie (d =1.32) sont plus lourds que l'eau (et done coulent), et l'if a une densite 
de 1. 

L'eau de mer et le lait ont une densite de 1.03 (on flotte micux dans l'eau de mer), le vin 
de 1.01 (on y flotterait legerement mieux aussi :-), et le petrole de 0.8 (c'est pourquoi on ne 
doit pas jeter de l'eau pour eteindre un feu d'hydrocarbure, car celui-ci flotte et la masse 
cnflammee se repand encore plus vite). 

Rcmarquons enfin que nous subissons egalement une poussee d'Archimede dans l'air, et la 
pression au niveau de nos jambes est effectivement plus elevee que celle au niveau de notre 
tete. Mais cette poussee est tres faible devant notre poids en raison de la faible masse volumique 
de l'air. 



46 



3.4 Forces de contact : frottement visqueux. 



3.3.8 Moment des forces de pression. 

Le moment d'une force de surface, notamment celle de pression, est obtenu en som- 
mant toutes les contributions de moments elementaires sur chaque element de sur- 
face. Ainsi le moment en un point A des forces de pression s'ecrira : 



M A {F P ) = J J AM A -pn dS (3.11) 



s 



Le calcul des moments est important par exemple pour les etudes de stabilite de bar- 
rage, de corps flottants, de bateaux .... On peut ainsi definir le centre de poussee 
C par : 
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CMA-pndS = (3.12) 
Par un raisonnement analogue au precedent, il est facile de montrer que : 



Le centre de poussee de la force d'Archimede exercee sur un corps est le centre 
de gravite du fluide deplace. 



Ainsi si le corps est homogene et totalement immerge, centre de poussee et centre 
de gravite du corps sont confondus. En revanche, si le corps a une partie emergee 
(notamment les bateaux) ce n'est plus le cas, et la stabilite du corps vis-a-vis de la 
rotation n'est plus assuree. 



3.4 Forces de contact : frottement visqueux. 



3.4.1 Mise en evidence : experience de Couette. 

Lorsque le fluide est en mouvement, on constate experiment alement qu'il apparait 
egalement une force tangentielle, en plus des forces de pression. L 'experience la 



plus simple permettant de le constater est l'ecoulement de Couette (Fig. 3.5) : un 
canal contenant du fluide initialement immobile est equipee d'une paroi superieure 
mobile. Lorsque Ton anime cette derniere d'un mouvement horizontal uniforme a la 
vitesse U on constate que : 



- le fluide est entraine par la plaque, avec une vitesse horizontale 

- la plaque subit une force tangentielle F opposee au mouvement. 



3.4 Forces de contact : frottement visqueux. 
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Fig. 3.5 - Experience de Couette : la plaque superieure entraine le fluide a cause 
des frottements visqueux. Le profil de vitesse est lineaire en regime permanent. 

Au bout d'un temps suffisamment long, la vitesse du fluide devient constante au 
cours du temps (on est en regime permanent) et on constate qu'elle varie lineairement 
de a Uq dans l'epaisseur du fluide. 

On peut conclure de cette experience que : 

- la plaque "adhere au fluide" puisqu'elle l'entraine dans son mouvement. 

- les couches horizontales de fluide frottent l'une sur l'autre puisque le mouvement 
se transmet de la plaque vers les couches plus basses. Cela signifie que chaque 
"lame" horizontale subit une force tangentielle de la part de ses voisines. Cette 
force tangentielle est appelee cisaillement [^} 

Autrement dit cette experience simple ^jdemontre l'existence d'une force tangentielle 
dans un fluide en mouvement causee par le frottement des couches fluides entre elles. 
On montre que la force par unite de surface est proportionnelle au gradient de 
vitesse dans la hauteur de l'ecoulement : 



F v dv x 

b dz 



3.4.2 La viscosite. 

Le coefficient de proportionnalite /1 est un parametre physique du fluide appele 
viscosite dynamique. Son unite S.I. est le Pa.s ou bien le kg.m~ :L s~ 1 , appele aussi 
Poiseuille^] et note PI. Son unite dans le systeme CGS s'appelle le Poise et est 



3 Le cisaillement existe aussi dans les solides, par exemple lorsque Ton deplace une surface 
du solide parallclcmcnt a cllc-meme. A l'inverse d'un fluide, qui coule lorsqu'il est soumis a du 
cisaillement externe, un solide tcnte de rcprendre sa forme initiale par une force elastiquc. 

4 On peut la realiser grossierement dans une baignoire en posant la main a plat a la surface de 
l'eau et en decrivant un mouvement circulairc. On sent au bout d'un certain temps que le fluide 
en profondeur s'anime lui aussi d'un mouvement de rotation 

5 Jcan-Louis Marie Poiseuille (1797-1869) : d'abord eleve de polytechnique, il etudie ensuite la 
medccine. Ses recherches mclangcnt lois physiques et physiologic, et il etudie notammcnt l'hcmo- 
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note Po. On verifiera que 1 PI = 10 Po. La viscosite de l'eau a temperature et 
pression ambiante vaut 10~ 3 PI, et celle de l'air 1.85 10~ 5 PI. La viscosite de l'eau 
vaut done un CentiPoise (cPo), et e'est pourquoi cette derniere unite est tres 
poupulaire. 

La viscosite varie notablement avec la temperature. Elle augmente avec T dans le 
cas d'un gaz et diminue dans le cas d'un liquide (vous pourrez verifier en regardant 
comment coule l'huile dans une poele chaude par rapport a une froide). La viscosite 
des liquides augmente egalement avec la pression, celle des gaz y est peu sensible. 



3.4.3 Origine microscopique. 

Dans l'experience precedente, nous voyons des filets de fluides qui se deplacent ho- 
rizontalement, aucun mouvement vertical n'est visible. Cependant a l'echelle micro- 
scopique les molecules se deplacent en tous sens, en plus du mouvement d'ensemble 



(voir Fig. 3.1). Sur la figure 3.5, une partie des molecules "hautes" et done rapides 
vient faire un tour vers le bas. Ce faisant elles choquent les molecules plus lentes et 
leur transfere de la quantite de mouvement. Reciproquement les molecules basses et 
lentes vont faire des incursions vers le haut et recuperent de la quantite de mouve- 
ment en choquant des molecules plus rapides. C'est ainsi que le deplacement de la 
paroi superieure finit par etre "ressenti" par les molecules du bas : de la quantite de 
mouvement est transferee dans l'epaisseur du fluide par l'agitation microscopique. 
Notons que ce mouvement vertical des molecules est a moyenne nulle, et done invi- 
sible a l'echelle macroscopique. 



3.4.4 Le nombre de Reynolds. 



Nous verrons ulterieurement que la prise en compte des forces visqueuses introduit 
une complexite supplementaire dans les equations du mouvement du fluide. C'est 
pourquoi il est pertinent de se demander dans quelle mesure on peut les negliger, 
en ne considerant que les forces de pression. L'experience montre par exemple que 
pour un ecoulement externe (fluide s'ecoulant autour d'un obstacle) les effets de 
la viscosite sont concentres dans une fine couche pres des parois solides, appelee 
couche limite^} Ceci est vrai pourvu que l'ecoulement ait une vitesse suffisante et 
que la viscosite soit sumsamment faible. 

En fait, l'epaisseur de la couche limite depend d'un nombre adimensionnel, appele 

dynamique (la circulation sanguine). Dans la lignee de ces travaux, il etudie le mouvement d'un 
fluide dans des tubes capillaires, et ecrit la fameuse loi de Hagen-Poiseuillc (Hagcn decouvrit la 
meme loi independamment), formulee plus rigoureusement un peu plus tard. Cette loi indique que 
le debit dans un tube est proportionnel a son diametre, a la difference de pression entre l'entree et 
la sortie, et inversement proportionnel a la viscosite du fluide (voir aussi section sur les pertes de 
charge) . 

6 ainsi que, dans une moindre mesure, dans le sillage qui est la zone d'ccoulcmcnt prolongcant 
la couche limite en aval de l'obstaclc 
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nombre de Reynolds^} fondamental en mecanique des fluides 



Le nombre de Reynolds est defini par : 




Re = PLV 


(3.13) 


C'est un donnee fondamentale d'un ecoulement. 





Dans cette definition, L est une longueur caracteristique de l'ecoulement (taille d'un 
obstacle, d'un recipient, diametre d'un tuyau, epaisseur d'un couche de fluide), v 
est une vitesse caracteristique de l'ecoulement. Le Reynolds Re peut etre interprete 
(entre autres) comme le rapport de l'energie cinetique et du travail des forces vis- 
queuses. 



Pour illustrer l'importance de ce nombre, la figure 3.6 presente l'ecoulement au- 
tour d'un cylindre pour 4 Reynolds differents (1.5, 26, 200 et 8000, respectivement, 
de gauche a droite et de haut en bas). On voit que l'ecoulement, au depart sy- 
metrique (Re=1.5), presente ensuite une dissymetrie amont-aval (Re = 26) et des 
re-circulations. Pour Re = 200, l'ecoulement se destabilise et devient instationnaire 
(au depart periodique), formant des tourbillons se detachant du cylindre alternati- 
vement de chaque cote|^] Enfin pour Re = 8000, l'ecoulement dans le sillage presente 
des mouvements erratiques, c'est l'apparition de la turbulence. 

De fagon generale, lorsque le Reynolds est tres petit devant 1, les effets visqueux 
sont importants dans tout l'ecoulement, c'est le cas par exemple des ecoulements 
ou bien tres visqueux (huiles, goudron, pates, plastique, miel), a tres faible vitesse 
(glaciers, manteau terrestre) ou d'objets microscopiques (bacteries, suspensions de 
particules). Dans ce cas les forces visqueuses jouent un role decisif dans l'ecoulement 
et ne peuvent etre negligees. 

A l'inverse lorsque le nombre de Reynolds est tres grand devant 1, et dans le cas 
d'un ecoulement externe (autour d'un objet, d'un vehicule. . .) les forces visqueuses 
seront concentrees au voisinage des parois solides et pourront etre negligees hors 
de cette zone. Dans le cas des ecoulements internes, cela est moins clair comme 
en temoigne clairemnt l'experience de Couette : nous avons vu que le frottement 



7 Osborne Reynolds (1842-1912) : sa contribution majeure fut en hydrodynamique a partir de 
1873. II etablit par une experience celebre les conditions de transitions d'un ecoulement laminaire 
vers un ecoulement turbulent. En 1886, il propose une theorie de la lubrification (ecoulements du 
fluide visqueux dans un espace de petites dimensions dans la direction perpendiculaire a l'ecou- 
lement). Trois ans plus tard, il produit le premier modele mathematique de la turbulence, dont 
les grands principes ont peu varie depuis (en particulier le fameux tenseur de Reynolds). II fut le 
second enseignant en ingenicric d'Anglcterre et possedait de fermes convictions quant a l'ensei- 
gncment de cette discipline : tout ingenieur, quel que soit son domainc, devait scion lui avoir une 
bonne connaissance des mathematiques, de la physique et de la mecanique classique. . . 

8 la photo est prise de plus loin pour visualiscr l'extension spatiale des tourbillons. Cette structure 
est appelee "allee de tourbillons de Von Karman". 
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Fig. 3.6 - Ecoulement autour d'un cylindre pour des Reynolds de 1.5, 26, 200 et 



8000, de gauche a droite et de haut en bas (d'apres Guyon et al. , 2001 ) 



visqueux etait dans ce cas le moteur du mouvement, et il serait paradoxal de negliger 
ces frottements, aussi faibles fussent-il ! 

Pour des ecoulements en conduite, le frottement sur les parois du tuyau conduisent 
a une perte d'energie mecanique du fluide, ce qui se traduit, pour une tuyauterie de 
section constante, a une chute de pression dans le sens de l'ecoulement. De meme, les 
singularites de la conduite, comme par exemple, un changement brusque de diametre, 
un coude, une fente. . . favorisent localement la dissipation d'energie par dissipation 
visqueuse et conduite a une perte d'energie mecanique. Cette perte est appelee perte 



de charge, et sera abordee dans la section |5.4.2| Comme nous le verrons, le modele 
du fluide parfait fournit malgre tout un resultat remarquablement simple, a travers 
la celebre formule de Bernoulli, que Ton peut corriger a posteriori pour prendre en 
compte les pertes de charge. 



3.4.5 Le modele de fluide parfait. 

Le modele de fluide parfait consiste a negliger Faction des forces visqueuses 
dans le fluide. 

II s'agit d'un modele. En realite, tous les fluides sont visqueux 
II sous-entend un mouvement du fluide non dissipatif. 



a Sauf l'helium superfluide, peii courant en genie des precedes 
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Nous verrons en effet que le modele de fluide parfait est associe a la conservation 
de l'energie totale (cinetique + potentielle) a la maniere d'un bille roulant sans 
frottement sur une glissiere. Un fluide parfait, s'il ne subissait aucune force exte- 
rieure, coulerait indefiniment une fois mis en mouvement. Ainsi, ce modele predit 
les comportements suivants : 

- un bateau ayant regu une impulsion de demarrage pourrait traverser l'ocean sans 
moteur 0; 

- un planeur planerait indefiniment meme sans courant ascendant ; 

- un fluide une fois mis en mouvement dans un tube horizontal pourrait s'ecouler 
sur des milliers de kilometres ; 

- les vagues, le son ne s'amortiraient pas. 

II convient de garder ces exemples a l'esprit lorsque Ton veut faire usage du modele 
de fluide parfait, mais il ne faut pas non plus penser que le fluide parfait est un 
concept inutile. Nous verrons dans la premiere partie de ce cours un grand nombre 
d'applications de ce modele, rendues possibles par la simplification drastique des 
equations du mouvement. Citons un exemple historique : le modele de fluide parfait 
a permis notamment de calculer de fagon analytique les forces de portances des ailes 
d' avion. 

Le probleme est qu'il est difficile d'etablir une regie autorisant ou non son emploi, 
tout simplement parce que pour un meme ecoulement, il se justifie ou non selon 
l'information que l'on cherche a calculer. Ainsi le modele fluide parfait fournit la 
portance d'une aile d'avion mais pas sa trainee. II fournit la vitesse de propagation 
des vagues ou du son, mais pas leur amortissement. On peut done enoncer les regies 
suivantes : 



Le modele de fluide parfait ne peut en aucun cas decrire : 

- les mouvements a faible nombre de Reynolds 

- les phenomenes dissipatifs, de perte, d'amortissement dans les ecoulements. 



Enfin, tres paradoxalement, il ne peut predire le phenomene de turbulence, qui 
pourtant apparait a tres grand nombre de Reynolds ! Plus generalement, la majorite 



mais pas tous ) des mouvements de fluides tourbillonnaires est provoquee par le 



frottement visqueux. 



9 Ce n'est pas tout a fait vrai meme si on neglige les forces visqueuses car le bateau perd de 
l'energie par les vagues qu'il provoque et qui se propagent loin de lui 

10 Les mouvements tourbillonnaires de l'atmosphere sont causes par la force de Coriolis. Des 
variations de concentrations dans des melanges peut egalement engendrer des tourbillons. 
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Fig. 3.7 - Contraintes normale et tangentielle. 



3.5 Ecriture tensorielle des forces de contact 



Nous avons vu que tout volume de fluide V est soumis a des forces de contact de 
la part du fluide environnant, appliquee sur sa frontiere S. Sur chaque element de 
la surface dS (Fig. 3.7[ ), la force exercee peut etre decomposed en une composante 
normale (de pression) et une composante tangentielle (dite de cisaillement), toutes 
deux proportionnelles a l'element de surface : 



dF = cr n dS + cr t dS 



(3.14) 



Le vecteur cr = dF/dS est appelee contrainte (en N/m 2 ). 

Pour decrire la contrainte en un point de la surface S dans un repere (x, y, z), il faut 
connaitre ses trois composantes pour toutes les orientations possible du vecteur n. 
En fait on montre qu'il suffit de connaitre cr seulement pour les 3 orientations de n 
selon x, y ou z, ce qui fait 9 composantes, que Ton peut organiser en une matrice 
3x3, appelee tenseur des contraintes, et notee : 



a 



j-y 



&yx ®yy 



a 



-■-il 



xz 



Ozz 



(3.15) 



n En fait, la compressibilite du fluide engendre aussi une force normale d'origine visqueuse. 
L 'etude detaillee des fluides reels compressibles n'entrant pas dans le cadre de ce cours, nous 
pourrons retenir contrainte normale = pression, contrainte tangentielle = visqueuse. 
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Fig. 3.8 - Illustration des composantes du tenseur des contraintes. 



La premiere colonne est la contrainte lorsque n est suivant Ox, la seconde lorsque n 
est suivant Oy, etc. . . Ou encore, a yx est la composante suivant Oy de la contrainte 



sur une surface dont la normale est suivant Ox (cf. Fig. 3.8) 



On notera que a xx , a yy , a zz representent les contraintes normales et les termes non 
diagonaux des contraintes tangent ielles. Enfin, et c'est l'interet de cette formulation, 
on peut montrer que la contrainte sur un element de surface de normale n peut 
s'ecrire simplement : 



dF = 

CT= ds =(Tn 

C'est un produit matrice-vecteur. Dans le cas du fluide, on a pour usage de "sortir" 
la pression du tenseur des contraintes, et de definir ainsi un tenseur des contraintes 
cf v uniquement lie a la presence des forces visqueuses. On aura done : 



cr = — = -pn + cr v .n (3.16) 
En integrant sur toute la surface S, nous pouvons done enoncer : 
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La force de contact exercee sur un volume de fluide V delimite par une surface 
S s'ecrit : 

F = J J -pn dS + J J W v .n dS (3.17) 

s s 



v v 

Forces de Forces visqueuses 

pression 



Dans le cadre du modele dit "de fluide parfait", on neglige la force vis- 
queuse (<r v = 0). 



Chapitre 4 

Equations du mouvement d'un 
fluide 



Dans le chapitre, precedent, nous avons dresse la liste des forces agissant sur un 
fluide, et il nous est maintenant possible de preciser la forme des equations de bilan 



de quantite de mouvement (2.14), et d'energie (2.17) en remplagant F cxt par le 



poids (3.1) + la force de contact (3.17). Nous differons volontairement l'expression 



de la contribution visqueuse, d'une part pour simplifier l'expose, et d'autre part pour 
obtenir un jeu d'equations independants du modele de fluide utilise. 

Nous presenterons dans un premier temps les equations de conservation sous forme 



de bilan volumique, immediatement deductible du chapitre |2.4[ Ensuite nous pre- 
senterons les formes locales de ces equations, sous formes d'equations aux derivees 
partielles, traditionnellement exposees dans la plupart des textes de mecanique des 
fluides. Enfin nous traiterons le cas particulier, mais fondamental des machines 
fluides, en relation avec le cours de thermodynamique des systemes ouverts. 



4.1 Sous forme de bilans volumiques 



4.1.1 Conservation de la masse. 



L'equation de conservation de la masse ne contient pas de termes source, et a deja 



ete presentee a la section 2.4.1 (Eq. 2.13 p. 31) 



II est interssant de rappeler sa forme, ou bien en regime permanent (d/dt = 0)) mais 
aussi pour les fluides incompressibles (p = C te ) sur lesquels nous reviendrons au cha- 
pitre ([5]). Nous appelons S e et S s les surfaces par lesquelles le fluide respectivement 
entre ou sort : 
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4.1 Sous forme de bilans volumiques 



Conservation de la masse en regime permanent d'un fluide dans un volume V 
fixe comportant une surface d'entree S e et une surface de sortie S s : 

-pv.n dS = / / pv.n dS (4.1) 



Masse entrante Masse sortante 



4.1.2 Conservation de la quant it e de mouvement. 

Nous conservons pour l'instant les forces de frottement visqueux, 1' approximation 
fluide parfait faisant l'objet du chapitre suivant. 

La forme generate de l'equation de conservation de la quantite de mouvement, pour 
un volume de fluide V fixe dans l'espace, s'ecrit : 



dt 




pvdV + / / (pv)v.n dS 



V s 
' ^ — 



Variation de qdm qdm sortante - 
dans V qdm rentrante 



— pn dS + 




pg dV 



v 



cr v .n dS 



(4.2) 



v 

Forces de 
pression 



Poids 



Forces visqueuses 



Ici encore il est interessant de considerer le cas du regime permanent avec une surface 
d'entree S e et une surface de sortie S s . On voit clairement que l'integrale de flux 
convectif est limitee aux surfaces d'entree et de sortie. Par ailleurs, l'integrals des 
forces de pression s'etend sur une surface fermee et en vertu du theoreme de la 
normale, on peut lui retrancher la constante p a tm- On a done : 



J J (pv)v.n dS = J J -(p - Patm)n dS + J J cr v .n dS + Mg 



Se~\~Ss 



L'interet de soustraire p atm vient de ce que la pression sur certaines surfaces (laterales 
ou en sortie) est bien souvent la pression atmospherique, si bien que les integrates 
de pression sur ces surfaces disparaissent de la formulation. La formule ci-dessus est 
utile en soi dans certains problemes (notamment les ecoulements fluviaux). 



4.1 Sous forme de bilans volumiques 
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Elle a un autre interet tres important : le calcul des forces exercees par le ffuide 
sur une structure solide avec laquelle elle est en contact. On considere la situation 



generale representee figure |4.1| : un volume de fluide est compris entre une surface 
d'entree S e , une surface de sortie S s , et les parois d'un solide, dont les autres faces 
sont soumises a la pression atmospherique. On cherche la force exercee par le fluide 
sur le solide. En vertu de la loi de Faction et de la reaction, elle s'ecrit : 



^■solide 



seche 




^■solide 



seche 



Fig. 4.1 - Fluide en contact avec un solide. 



FFluide/Solide — — F So lide/Fluide 



—pn dS 



cr v .n dS 



• 2-mouillcc 



^mouillc 



(4.4) 



En generale on s'interesse a la force totale exercee sur le solide, non seulement celle 
exercee par le fluide, mais aussi celle exercee par l'atmosphere sur ses parois seches 
SUche- Cette derniere s'ecrit : 



Air/Solide 



Patm n solide dS 



(4.5) 



Mais en appliquant le theoreme de la normale au solide, on peut ecrire : 



"Patmllsolidc 

dS 



Patmllsolidc dS 



-p atm n dS 



•Sscchc 



(4.6) 
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4.1 Sous forme de bilans volumiques 



ou la deuxieme egalite resulte de n so iid e = — n sur .Smouiiiee- Physiquement, la relation 
ci-dessus signifie que le solide repercute sur le fluide la force de pression qui lui est 



appliquee par l'air. On en deduit done la force totale, en sommant (4.4) et (4.5) : 



Totale/Solidc 



-(p-Patm)n dS + J J cr v .ndS 

' imouillcc 'S'mouillee 



(4.7) 



On utilise maintenant le bilan de quantite de mouvement (4.8), en decomposant les 



integrates de pression et de frottement sur S e + S s et S^uiiiee : 

J J (pv)v.n dS = J J -(p- Patm)n dS + J J a v .n dS 



S e -\-S s 



+ 

"Smouillcc 



Se~\~Ss 

-(p-Patm)n dS 



^ v .n dS+Mg 



(4- 



"^mouillc 



Totale /Solid c 



d'apres (4.7) 



On en deduit done le resultat fondamental : 



La force totale exercee sur un solide en contact avec le fluide d'un cote, et 
ratmosphere de l'autre s'ecrit 

F T otaic/Soiidc = J J ~{p- £>atm)n dS + J J W v .n dS - J J (pv)v.n dS + Mg 

(4.9) 

En general la contribution visqueuse a l'entree et a la sortie est negligeable. 



Cette formule est tres pratique car elle ne fait intervenir que les surfaces d'entree 
et de sortie, qui sont en general geometriquement simples, de teles sorte que les 
integrates sont faciles a calculer. 



4.1.3 Conservation de l'energie. 



Pour etablir l'equation de conservation de l'energie d'apres (2.17), il nous faut 
d'abord calculer la puissance des forces exterieures W. La puissance est le pro- 
duit scalaire de la force par la vitesse et il suffit done d'integrer toutes les puissances 
infinitesimales, sur le volume V pour les forces de volume et sur la surface S pour 
les forces de contact, d'ou : 



4.1 Sous forme de bilans volumiques 
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W= JJJ pg.vdV + J J -pn.vdS + J J v.(a v .n) dS (4.10) 

V s s 

V v ' V v ' V v ' 

Puissance du Puissance des Puissance des 

poids. forces de forces de 

pression. frottement 
visqueux. 



La conservation de l'energie d'un volume de fluide V fixe dans l'espace s'ecrit : 



d 
dl 




p[u + 



dV + 



v 



Variation 
(E interne + 
E cinetique) 




pg.v dV + 



p[u + 



v.n dS 



(E interne + 
E cinetique) 
sortante - 
rent rant e 

-pn.v dS + 



v.(ovn) dS + Q 



(4.11) 



Puissance du 
poids 



Puissance des Puissance des 

forces de pression forces visqueuses 



W puissance des 
forces exterieures 



4.1.4 Complement : theoreme de l'energie cinetique. 

II ne s'agit pas d'un principe de conservation en plus, mais d'une consequence de la 
conservation de la quantite de mouvement. Nous demontrons dans P annexe [Pique : 
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4.2 Conditions aux limites. 



V s 



V v ' ^ 



Variation E E cinetique 

cinetique sortante - 

rentrante 




pg.v dV +11 — pn.v dS + J J v.(cr„.n) dS 

. (4-12) 



v s 



W puissance des 
forces exterieures 




+ / / / pdivv dV - $„ 



v 



Wi n t puissance 
des forces 
interieures 



4.2 Conditions aux limites. 



La resolution des equations du mouvement necessitera la connaissance de la vitesse 
du fluide sur frontieres du domaine. Bien souvent ces limites sont des parois solides 
ou des interfaces avec d'autres fluides (liquides ou gaz). 

Nous ne traiterons que le cas le plus courant des parois solides, qui en aucun cas 
ne peut etre traverse par du fluide (a l'exception des parois poreuses dont nous ne 
parlerons pas ici) . Par consequent la composante normale de la vitesse du fluide sur 
une paroi solide est toujours egale a la vitesse normale de cette paroi : 



V'H-solidc VgQijdQ.ngojidg 



Cette condition est connue sous le nom de condition d'etancheite. 

Dans le cas d'une paroi fixe, cette condition devient : 



v - n solide — 



Concernant la composante tangentielle, comme en temoigne par exemple l'experience 
de Couette, la composante tangentielle du fluide est egale a la vitesse tangentielle 
de la paroi, et done nulle pour une paroi fixe [j] En combinant cette condition avec 
la condition d'etancheite ci-dessus, on a pour un fluide reel : 



1 On pourra s'en convaincre en passant sa main par la vitre d'une voiture et en l'approchant 
lentement de la portiere. On sentira une diminution progressive de la vitesse. 



4.3 Application aux ecoulements en tuyauterie. 
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V — V so ii de 

Cela est du au fait qu'une paroi solide est constitute au niveau microscopique de 
petites asperites qui interdisent un mouvement d'ensemble des molecules. C'est pour 
cette meme raison que le solide subit une force tangentielle de la part du fluide. Aussi 
cette condition est-elle liee au caractere visqueux du fluide. 

Par contre si Ton neglige ces frottements visqueux, il faut etre coherent et laisser le 
droit au fluide de glisser librement sur la paroi ! Aussi dans le cadre du modele de 
fluide parfait, seule la condition d'etancheite sera ecrite. 

On retiendra done : 



La condition sur la vitesse au niveau 


d'une paroi solide est : 




- pour le modele de fluide parfait, 






V.n so iide 


VsolidcHsolide 


(4-13) 


- pour le modele de fluide reel (avec 


frottements visqueux) 




V 


^solide 


(4.14) 



Pour les ecoulements externes, on utilisera souvent des conditions rejetees a Finfini 
de l'obstacle du type : 

limv = V 00 (4.15) 

r— >oo 

ou r est la distance a l'obstacle. Voo est par exemple l'oppose de la vitesse d'un 
vehicule, lorsque Ton etudie l'ecoulement dans un referentiel lie a ce dernier. 

Pour les ecoulements fermes (en general en tuyauterie), on connait en general la 
vitesse a l'entree, et la pression a la sortie (par exemple pression atmospherique si 
l'ecoulement debouche a l'air libre). 



4.3 Application aux ecoulements en tuyauterie. 
4.3.1 Preliminaire. 

II est interessant d'un point de vue pratique, et pour faire le lien avec le cours 
de thermodynamique, de reprendre les equations de conservation sur un tube de 
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4.3 Application aux ecoulements en tuyauterie. 



courant, en faisant l'hypothese d'un ecoulement piston, afin de pouvoir supposer les 
grandeurs homogenes sur les sections d'entree et sortie. On suppose de plus les parois 



solides. Le systeme est represented figure 4.2), on retrouve la section d'entree S e , la 



section de sortie S s et la surface laterale est constitute par un solide en contact 



avec le fluide. Nous allons ecrire l'expression generate (2.12) rappelee ci-dessous 



respectivement pour la masse, la quantite de mouvement et l'energie. 



dG 
~dt 



g e v e S e - g s v s S s +R + (t) - R (t) 



Ge 



G, 



(4.16) 




Fig. 4.2 - Tube de courant. 



4.3.2 Conservation de la masse. 



n, n s 



II n'y a pas de terme source dans la conservation de la masse (R + = R = 0) et 



1' equation (2.12) s'ecrit done, avec G = M et g = p 



La conservation de la masse sur le tube de courant s'ecrit alors simplement : 

dM 



dt 



p e v e S e - p s v S S S 

Me M a 



(4.17) 



Dans le cas du regime permanent, le debit massique est conserve dans le 
tube de courant : 



M s = M, 
p e v e S e = p s v s S s 



(4.18) 
(4.19) 



Par unite de temps, il rentre autant de matiere dans le tube qu'il n'en sort, 
autrement dit le tube n'accumule pas de matiere. 



Ce resultat tres intuitif est l'un des fondements du genie des precedes. 



4.3 Application aux ecoulements en tuyauterie. 
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4.3.3 Conservation de la quantite de mouvement. 



En ecrivant (2.12) avec G — P, g — pv, R + — R = ^ F ext , nous obtenons 



dP 

— = p e v e S e v e - p s v s \ s S s + ^2 F cxt (4.20) 
= M e v e - M s w s + F ext (4.21) 

On notera que ces relations sont des relations vectorielles, ce qui est important 
notamment quand les surfaces d'entree et de sortie sont orientees differemment (voir 
exemple du tube en U ci-dessous). 

On remarquera que puisque F ex t est la somme des forces exercees par l'exterieur 
sur le volume de fluide, — F ext est la force exercee par le fluide sur l'exterieur, 
et apparait dont comme la difference des flux de quantite de mouvement entrant 
et sortant. Ce dernier resultat est d'une utilite tres pratique, par exemple pour 
determiner les efforts exercees par l'ecoulement du fluide sur une tuyauterie coudee^] 
ou a section variable (voir TD). II est egalement a la base du principe de propulsion 
(reacteur, helice . . .) ou Ton accelere le fluide qui rentre grace un apport d'energie 
(en brulant du combustible). 

Ces derniers points peuvent etre mis en evidence tres simplement en explicitant les 



integrales sur S e et S s dans l'expression de la force totale sur le solide (4.9) obtenue 
precedemment. En admettant que les forces visqueuses sur les surfaces d'entree 5* e 
et de sortie S s sont negligeables (elles le sont dans la pratique), on obtient 



La force totale exercee par un ecoulement permanent + l'air exterieur sur un 
tube solide s'ecrit 



F T otai e /Soiide = [M e v e + (p e - v^S^e ~ [M s v s + (jp s - p atm )S' s ]n s + Mg 

= {p e V 2 e +p e - p at m)S e Il e - (p s V 2 s +p s - p at m)S s Il s + Mg 

(4.22) 

ou n e et n s sont les normales aux surfaces d'entree et de sortie, orientees dans 
le sens de l'ecoulement. 



En l'absence de mouvement, la force exercee sur la structure est compensee par le 
poids du fluide + les forces de pressions a l'entree et a la sortie. En presence de 
mouvement, la structure subit une force additionnelle egale a la perte de quantite 
de mouvement du fluide entre l'entree et la sortie. 



2 On pourra en faire l'experience avec un tuyau flexible (de douche par exemple) qui a une 
facheuse tendance a se promener lorsquc Ton ouvre le robinct 
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4.3 Application aux ecoulements en tuyauterie. 



La formulation ci-dessus est particulierement interessante lorsque la tuyauterie de- 
bouche a l'atmosphere, auquel cas p s = p a tm- 



Application : considerons un tube en U de section S (Fig. 4.3) dans lequel un 
ecoulement d'air rentre a la pression atmospherique p atm et a la vitesse v = ve x . 
Nous supposerons le fluide incompressible de telle sorte que p e = p s = p. 



n 



n 




Fig. 4.3 - Force exercee par un ecoulement sur un tube en U 

Dans ces conditions, la conservation de la masse en regime permanent impose que 
la vitesse de sortie soit egale a — v e x . En negligeant le poids de Pair dans le tube, la 



formule (4.22) fournit : 



FTotale/Solidc — 2pV 2 Se x 

et est done orientee dans le sens de e x . 

On voit done que la force totale exercee sur la tuyauterie est independante de la 
pression atmospherique. Ce resultat est vrai pour tout ecoulement interne pour lequel 
la paroi du tube est soumise a la pression atmospherique. 



4.3.4 Conservation de l'energie. 



Nous allons ecrire (2.12) avec G = U + K, g = p(u + v 2 /2) et R + - R~ = W + Q. 



Le calcul de Q releve du cours de transfert thermique et sera laisse tel quel ici. Nous 
pouvons en revanche calculer W, puisque nous connaissons les forces exterieures : 
poids, force de pression et de frottement visqueux. Nous allons generaliser un peu le 
cas du tube de courant de la fagon suivante : 

- nous rajoutons une surface S u mobile a l'interieur du tube, correspondant par 
exemple aux pales d'une helice ou d'une turbine. 

- nous supposons que S^t est une paroi solide, qui en vertu de la condition d'adhe- 
rence, est bien une ligne de courant. 

Ces hypotheses minimales permettent de representer un grand nombre de compo- 



sants des reseaux de fluide. II est facilement verifiable que l'equation (2.12) reste 
done valable en presence de la surface mobile S u . 



4.3 Application aux ecoulements en tuyauterie. 
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Demonstration : l'ajout de S u qui est une surface mobile pose une legerc dimcultc supplcmen- 
taire car nous avons suppose jusqu'alors que le volume de bilan V etait un volume fixe, ce qui 
n'est plus le cas ici. Comme nous 1'avons mentionne page 24 le terme convectif J J gv.n dS 



doit etre generalise en J J g(v — w).n dS, ou w est la vitesse locale de la surface S. Ici le 

s 

probleme se pose seulement sur S Ul qui est la seule surface mobile. Mais cette surface est 
etanche et par consequent, nous avons v.n = w.n en chaque point de S u , et le terme convectif 
associe est done nul. 




Fig. 4.4 - Schema generique d'un machine fluide 



En revanche la presence de S u influe bien sur sur le calcul de W, puisque cette 
piece mobile est justement presente pour fournir (ou prelever) de l'energie au fluide. 
Nous admettrons que la puissance des forces de frottement visqueux sur les sections 
d'entree et de sortie sont negligeables. De plus la puissance des forces de contact 
sur Siat est nulle, puisque la vitesse y est nulle de par la condition d'adherence (cf. 



section 



4.2). Ensuite, nous notons W u la puissance des forces exercees par la 



pale sur le fluide et nous supposerons la puissance du poids nulle pour simplifier 
l'expose^} Dans ces conditions, la puissance des forces exterieures (4.10) s'ecrit : 



3 ce qui est vrai si les sections d'entree et de sortie sont de faibles dimensions, a la meme hauteur, 
et que l'extension verticale des pales est petite. On consultera l'annexe [K\ pour la prise en compte 
de ce terme. 
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4.3 Application aux ecoulements en tuyauterie. 



W = J J -pv.ndS + J J -pv.ndS + J J (W v .n).v dS 

Se Ss Se~\~Ss 



~*V — 



" p e v e S e = -p s v s S s ~ (admis) 

+ / / (— pn + a v .n).\ dS + (-pn + <f„.n).v dS 



Sla,t S u 



= car v = 
sur Si a t (condition 
d'adherence) 



W u par definition 
///(pg).vdy (4.23) 




v 



suppose nul ici, 
cf. annexe lAl 



soit 



W = p e v e S e - p s v s S s + W u (4.24) 

Les deux premiers termes sont les puissances des forces de pression exercees par le 
fluide a l'entree et a la sortie sur le fluide interne a la machine. II est important de 
bien comprendre que W u et W different l'un de l'autre par ces deux termes. 



Dans ces conditions, (2.12) s'ecrit, avec G = U + K , g = p{u + v 2 /2) : 



d (U + K)= Pe (u e + ^] v e S e ~Ps(u s + ^\ v s S s 



dV ' \ 2 J r \ 2 

+p e v e S e - p s v s S s + W U + Q 



Mais en rappelant que le debit massique s'ecrit M = pvS, le bilan se re-ecrit : 



d .(U + K) = M e (u e + ^ + V 4) -M s (u s + ^ + 4) + -W.. + Q 



dt V Pe 2 J \ p s 2 

ce qui fait apparaitre l'enthalpie massique h = u + pj p, d'oii finalement 



^-(U + K) = M e (h e + V 1\ - M s (h s + V l) +W U + Q (4.25) 



Ici nous n'avons pas pris en compte le travail des forces de pesanteur. Dans le cas le 
plus general oil la pesanteur joue un role notable, notamment si les sections d'entree 



4.3 Application aux ecoulements en tuyauterie. 
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et de sortie sont a des altitudes tres differentes, nous montrons en annexe [A] (la 
demonstration est sans dimculte particuliere mais un peu longue) que ce resultat se 
generalise en : 



Le premier principe en systeme ouvert s'ecrit 



j t (U + K + E p ) = M e (h e + V i+ gz)j -M s (h s + V ^+ gz)j +W U + Q (4.26) 



ou E p est l'energie potentielle de pesanteur du fluide dans la machine. 
Insistons sur les points suivants : 

- l'enthalpie provient des forces de pression a l'entree et a la sortie, que nous avons 
separe volontairement du travail utile, tout simplement parce que seul ce dernier 
est interessant dans ce genre de probleme. 

- le terme instationnaire dll/dt est une variation d'energie interne et non pas 
d'enthalpie comme parfois ecrit a tort. A ce titre, il est done injustifie de parler 
de "bilan enthalpique", comme e'est malheureusement souvent le cas en genie des 
precedes. 

- la puissance W u intervenant ici, dite puissance utile, est la puissance des seules 
forces exercees par les pales a l'interieur de la machine, car les puissances des 
forces de pression a l'entree et a la sortie sont deja integrees dans les enthalpies. W u 
est nulle dans plusieurs cas pratiques (vanne de detente, echangeurs de chaleur), 
positive si le fluide recoit du travail de ces pales (compresseur, pompe), negative 
si le fluide fournit du travail (turbine, eolienne, moulins). 



4.3.5 Theoreme de l'energie cinetique. 



II est egalement interessant d'ecrire le theoreme de l'energie cinetique (4.12) dans 
cette configuration. La calcul est assez simple, car nous venons de calculer la puis- 
sance des forces exterieures. Nous obtenons : 
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4.3 Application aux ecoulements en tuyauterie. 



Le theoreme de l'energie cinetique pour une machine fluide s'ecrit, dans le cas 
general : 



*. {K + E r ) = M e gj + f + 92e ) - M, + f + 9 z. ) + W. 




+ pdivv - (4.27) 



v 





Puissance des 
forces interieures 



Nous voyons que cette equation ressemble beaucoup au premier principe, mais il ne 
faut cependant pas les confondre. La derniere equation fait intervenir la puissance 
des forces interieures, absente dans le premier principe. 

- Le premier terme JJJ pdivv dV est la puissance des forces de pression interne 

v 

(nulle en incompressible) et correspond a de l'energie elastique emmagasinee par 
le fluide dans la machine. II peut etre positif (cas d'une detente) ou negatif (pour 
une compression). 

- Le terme en revanche est toujours positif et represente la puissance des forces 
de frottement visqueux a l'interieur du fluide. II s'agit d'une dissipation irre- 
versible d'energie. Nous verrons ulterieurement que ce terme est fondamen- 
tal notamment pour decrire les pertes d'energie dans les ecoulements en 
conduite. 



II est enfin interessant de soustraire cette equation au premier principe (4.26) pour 
obtenir les variations d'energie interne du fluide : 



= M e u e - M s u s - JJJ pdiv v dV + + Q (4.28) 
v 

Nous voyons que la dissipation irreversible d'energie toujours positive, peut par 
exemple : 

- augmenter l'energie interne du fluide U dans la machine (done sa temperature) 

- etre compensee par une cession de chaleur a l'exterieur {Q < 0) 

- augmenter l'energie interne u s du fluide sortant (done sa temperature) 

Ainsi le simple fait de mettre en mouvement un fluide dans un recipient ferme 
(M = 0) et calorifuge (Q = 0) echauffe ce fluide (pensez aux lubrifiants dans une 
voiture par exemple). 



On voit egalement que la puissance des forces interieures de pression a le meme 
effet dans le cas d'une compression, et l'effet inverse dans le cas d'une detente. Par 



4.4 Equations locales. 
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exemple dans le detendeur d'un refrigerateur, qui fonctionne en regime permanent, 
ce terme est negatif et est partiellement compense par un Q > (le detendeur regoit 
de la chaleur du compartiment du refrigerateur . . . et refroidit done les aliments qui 
y sont stockes). 



4.4 Equations locales. 

4.4.1 Contexte 

Jusqu'a present, nous avons ecrit des equations de bilan sur un volume V macro- 
scopique, fixe dans l'espace. II est possible d'en deduire des equations aux derivees 
partielles, dites locales, reliant les grandeurs de l'ecoulement en chaque pointj^] 

Leur interpretation directe pose souvent des problemes au debutant, et e'est pour- 
quoi nous avons reporte leur ecriture le plus loin possible. Nous verrons que sous 
certaines conditions, il existe une forme integree (sans aucune derivee ni integrate !) 
de ces equations, dite formule de Bernoulli, tres pratique dans un certain nombre de 
cas. 

Aussi, les equations locales seront peu utilisees dans cette premiere partie du cours, 
sauf : 

- pour demontrer le theoreme de l'energie cinetique enonce plus haut 

- pour demontrer la formule de Bernoulli presentee au chapitre suivant pour un 
fluide parfait. 

II faut bien comprendre cependant que la majorite des problemes de mecanique des 
ffuides ne peut etre resolu que par les equations locales et elles constituent done un 
outil indispensable dans le cas general. 

4.4.2 Obtention 

Leur ecriture se deduit des equations precedentes a l'aide des formules d'Ostro- 
gradski (cf. annexe [B| . Nous detaillons seulement l'equation de conservation de la 



masse : reprenons (|4.1|) : l'integrale de surface se transforme par (|B.3|) et la derivee 
a_ 
dt 



d commute avec l'integrale triple car V est un volume fixe. On obtient done : 




° P dV+ I I I div (pv) dV = 



dt 

V V 




Les cours de mecanique des fluides commencent souvent par developper ces equations. Nous 
preferons la demarche inverse dans la mesure ou les bilans sous forme globale nous paraissent plus 
intuitifs. 
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4.4 Equations locales. 



Cette equation est vraie pour tout volume V aussi petit fut-il et par consequent 
l'integrande est nul, so it : 



^ + div (pv) = 


(4.29) 


conservation de la masse sous forme locale. 





De la meme fagon, en utilisant (B.4) on demontre l'equation de conservation de la 
quantite de mouvement, dite "sous forme conservative", car elle decoule directement 
du bilan ecrit plus haut : 

<9(pv) — 

— h div (pv <g> v) = - gradp + pg + div cr v (4.30) 

at 

On remarquera l'apparition de deux objets mathematiques inusuels 

- le produit v ® v est le tenseur de terme principal V(Vj (ecrivez la matrice a titre 
d'exercice) 

- l'operateur div (qui fournit un vecteur, a ne pas confondre avec div d'un vecteur 
qui fournit un scalaire) s'applique a un tenseur et donne un vecteur dont chaque 
composante est la divergence de chaque ligne de la matrice. 

On utilise plus souvent la forme non conservative de cette equation. On montre que 
le premier membre peut etre transforme, en utilisant la conservation de la masse, 
par : 



<9v — 
p-7^7 + p(v. V) v = - grad p + pg + div cr v (4-31) 

conservation de la quantite de mouvement sous forme locale. 



Le produit (v. V) est formellement le produit scalaire du vecteur vitesse par le vecteur 
nabla et est done (en coordonnees cartesiennes) l'operateur v x -^ + v y -^ + v z -j^. Par 
consequent le terme (v.V)v est cet operateur applique a chaque composante de v, 
soit : 



+ v. 



dv x 



+ v z 



J X r% I U V r\ i "Z r\ 

ox oy oz 



(v.V)v 



+ v v 



dv v 



+ v. 



ox oy oz 



dv. 



dv z 



dv z 



v x ^— +v v — 1 + v z — 
ox ' oy oz 



4.4 Equations locales. 
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29), (4.31) sans les forces de frottement visqueux forme les 



Le jeu d'equations 
"equations d'Euler 



Enfin l'equation de conservation de l'energie peut etre obtenue a partir du premier 



principe (4.11), puis par utilisation de la formule (B.3). On obtient, sous forme 



conservative : 



d[p(u + v 2 /2)} 
dt 



+ div [p(u + v 2 /2)} = - div (pv) + div (W v .v) + pg.v - div q (4.32) 



ou q est une densite de flux de chaleur sortant (Q = — J J q.n dS), dont la specifica- 

s 

tion est l'objet du cours de transfert thermique. Ici encore le membre de gauche 
peut etre transforms? pour obtenir : 



Q( u _|_ v 2 _ 
p — h pv. grad [u + v 2 /2) = — div (pv) + div (<f„.v) + pg.v — div q 

(4.33) 

expression du premier principe sous forme locale. 



L'equation de l'energie peut etre ecrite sous des formes diverses, Notamment en lui 
soustrayant, comme nous l'avons fait sous forme globale, le theoreme de l'energie 



cinetique. Le lecteur interesse pourra se reporter au livre de Bird et al. (1960) 
page 322. 



4.4.3 Un jeu d'equations complet ? 



A ce stade il est pertinent de se demander si le jeu d'equations (4.29), (4.31), (4.33) 
est sufnsant pour resoudre un probleme pratique. Pour cela comptons les equations 
et les inconnues : 

- nous avons 5 equations : 1 pour la masse, 3 pour la quantite de mouvement car 
c'est une equation vectorielle, 1 pour l'energie. 



5 Leonhard Euler (1707-1783), mathematicien Suisse : initialement destine par sa famille a une 
carriere ecclesiastique, il se consacrera avec succes aux sciences. II en effet l'auteur de nombreux ou- 
vrages en mathematiques, astronomie, mecanique, physique, optique ou encore acoustique. II sera 
surtout connu pour ses travaux mathematiques, dans le domaine de l'analyse (etude de fonctions 
et calcul infinitesimal). II fut rclcve de Jean Bernoulli et ami de son fils Daniel, egalement a l'ori- 
ginc d'unc formule cclcbrc en hydrodynamique. II inspirera de nombreux chercheurs, notamment 
Lagrange. 

6 On admettra pour l'instant que le tcnscur des contraintes visqucuses 7f v et la densite de flux 
de chaleur q n'introduisent pas de variables supplementaires. 
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4.4 Equations locales. 



- nous avons 6 inconnues : la masse volumique p, les 3 composantes de la vitesse v, 
la pression p et l'energie interne u. 

II manque done une equation ! Cette equation manquante est l'equation d'etat 
reliant par exemple u, p et p. 



Les equations de conservation doivent etre completees par une equation d'etat 
du fluide pour que le probleme puisse etre resolu. 



4.4.4 Cas du fluide incompressible 

II existe une simplification majeure si Ton suppose le fluide incompressible. Dans 
ce cas la masse volumique p est supposee constante au cours du temps et dans tout 
I'espace. Cette variable disparait done de l'ensemble des inconnues et mieux encore, 
on voit alors que l'equation de l'energie est decouplee des equations de conservation 
de la masse et de la quantite de mouvement. 



Un ecoulement de fluide incompressible est defini par : 

p = C tc dans le temps et I'espace 
qui fait office d'equation d'etat. 

II peut etre decrit uniquement par les equations de conservation de la masse et 
de la quantite de mouvement. 

L'equation de conservation de la masse se reduit alors a : 

divv = (4.34) 



Ce decouplage de l'equation de l'energie est en fait vrai dans une classe de problcmcs plus 
large, ou Ton peut trouver une relation entre p et p (on parle de fluides "barotropes"). Par 
exemple pour un gaz parfait en ecoulement isentropique, on pourra ecrire pj ' p 1 — constante. 

II reste une question majeure : quand peut-on dire qu'un fluide est incompressible ? 
II faut considerer la question avec la meme approche que pour la prise en compte ou 
non des forces visqueuses. Tous les fluides sont en fait compressibles ! 

Ici l'intuition peut etre trompeuse : on sent bien qu'un liquide par exemple est 
moins compressible qu'un gaz, et Ton pourrait etre tente d'appliquer l'approxima- 
tion incompressible aux liquides mais pas aux gaz. Ce serait une erreur comme en 
temoignent les exemples suivants : 



4.4 Equations locales. 
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- l'ecoulement d'air autour d'un avion volant a faible vitesse peut etre considere 
comme incompressible 

- l'etude des explosions sous-marines serait grandement faussee par l'hypothese in- 
compressible 

- le son se propagerait a une vitesse infinie dans un fluide incompressible, liquide 
ou gaz. 

De meme que la prise en compte des effets visqueux depend du nombre de Rey- 
nolds, la prise en compte de la compressibilite depend de la valeur d'un nombre 
adimensionnel : 

On definit le nombre de MachQpar : 

Ma = - (4.35) 

c 

ou v est un ordre de grandeur de la vitesse de l'ecoulement, et c la vitesse 
du son dans le fluide. En general, on pourra raisonnablement considerer qu'un 
ecoulement est incompressible pour : 

Ma < 0.3 



a Ernst Mach (1838-1916) : scientifique autrichien. Apotre du positivisme (la science ne 
soit se baser que sur l'observation experiment ale et la sensation des phenomenes), son nom 
reste celebre par son approche des ecoulements supersoniques (Ma > 1), dont il pressent les 
proprietes particulieres. Avant Einstein, il est le premier a rejetcr le caractere absolu du temps 
et de l'espace, initiant ainsi la theorie de la relativite gencralc. II fut egalement l'un des plus 
fervents opposants a Boltzmann, et plus generalement a l'existence de l'atome. 



II existe bien sur des exceptions a cette regie. Par exemple, si l'on veut etudier la 
propagation du son, on ne pourra pas negliger la compressibilite du fluide, bien que 
les vitesses mises en jeu soient infimes. 

Quelques ordres de grandeurs de la vitesse du son : 340 m/s dans l'air 1500 m/s 
dans l'eau. 



Chapitre 5 

Mouvement du fluide parfait 
incompressible 



5.1 Rappel des hypotheses et equations 

Comme nous l'avons vu aux chapitres precedent : 

- le fluide peut en general etre considere parfait (les forces de frottement visqueuses 
sont negligeables) a grand nombre de Reynolds et loin des parois solides. Au niveau 
des parois, le fluide doit etre libre de glisser (vitesse tangentielle non nulle). 

- le fluide peut etre considere incompressible pour des nombres de Mach faibles 
(pratiquement < 0.3). 

Dans ces conditions les equations du mouvement sont, sous forme globale : 



v.n dS = 



d 
di 



J pv dV + J J (pv)v.n dS = J J —pn dS 




(5.1) 

pg dV (5.2) 



v 



Sur un tube de courant, ces equations deviennent : 



v e S e = v s S s (5.3) 

dP 

— = ( yPe + pvl) - (p s + pvl) - F Fluide/So iidc + Mg (5.4) 

ou P est la quantite de mouvement totale du tube, M est la masse de fluide dans le 
tube, vS — V est le debit volumique en m 3 /s, et F F i uidc /soiidc est la force exercee 
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5.2 Formule de Bernoulli 



par le fluide sur les parois laterales du tube plus eventuellement les pales dans le cas 
d'une machine. On comprendra mieux desormais pourquoi les fabricants de pompes 
(destinees a mouvoir du liquide) parlent plutot de debit volumique que de debit 
massique. 

Enfin, les equations du mouvement locales s'ecrivent : 



divv = (5.5) 
p— + (v.V)v = -gradp + pg (5.6) 

On remarquera dans les trois jeux d'equations que le terme instationnaire (en d/dt) 
disparait de l'equation de conservation de la masse. En regardant par exemple le 
bilan sur un tube de courant, on voit qu'une difference de debit d'entree et de sortie 
ne saurait etre compensee par une accumulation de matiere dans le tube, puisque le 
fluide ne peut etre comprime. 

Ces equations comportent toujours des derivees partielles et sont a ce titre, "com- 
plexes". Une simplification remarquable, rendant 1'hydrodynamique accessible au 
plus grand nombre en a ete proposee Bernoulli] au 18 eme siecle. 



5.2 Formule de Bernoulli 

5.2.1 Hypotheses-Enonce 



Avec les hypotheses : 

- fluide parfait, 

- fluide incompressible, 

- regime permanent, 

/ v 2 \ 

la quantite I p + p h pgz est constante le long d'une ligne de courant. 



1 Daniel Bernoulli (1700-1782) : scientifique Suisse, appartint a une famille de scientifiques. Son 
pere Jean Bernoulli fut un mathematicien brillant (anecdotiquement on lui doit le symbole g 
pour la pesanteur), dont Euler fut l'eleve. Son oncle Jacques et son frere Nicolas furent egalement 
mathcmaticicns. Daniel Bernoulli publie ses travaux sur l'liydrodynamique et la theorie cinetique 
des gaz en 1738 (il donna une version partielle de l'equation de Van der Waals avec un siecle 
d'avance). II s'interesse egalement a la theorie de l'clasticitc et aux vibrations de cordes et de 
tuyaux avec Euler, dont il fut l'ami proche, et d'Alembert. II fut egalement professeur de botanique, 
d'anatomie et de philosophie ( ! !). 



5.2 Formule de Bernoulli 



77 



5.2.2 Demonstration 



La demonstration, qu'on pourra oublier, est presentee ci-dessous. 



En regime permanent, le terme ^ disparait de l'equation (5.6 1. Par ailleurs en utilisant la 
formule (B.l I de 1' annexe |b) et en divisant tout par p, on obtient : 

,,2 



grad 



rot v x v 



gradp 

P 



grad gz 



en remarquant de plus que, avec 1'axe Oz oriente vers le haut g = — grad gz. Effectuons le 
produit scalaire de cette equation par un element dM de ligne de courant : puisque les lignes 
de courant sont tangentes a v, les deux vecteurs dM et v sont paralleles, et le produit mixte 
(rotv x v).rfM est done nul. II rcstc done : 

grad — .dM + gradp .dM + grad gz.dM = 
2 p 

Comme le fluide est incompressible, p peut etre entre dans gradp et nous obtenons : 

,,2 



grad 



P 
P 



gz .dM = 



ce qui signifie que l'interieur du gradient nc varie pas le long de la ligne de courant. CQFD. 



5.2.3 Comment aires 



Cette formule est certainement la plus connue de la mecanique des fluides, la plus 
pratique, et aussi la plus dangereuse. Elle est pratique car on voit qu'elle ne presente 
plus aucune derivee partielle, contrairement aux equations locales. Elle est dange- 
reuse car comme nous le verrons, elle ne s'applique pas dans certaines configurations 
ou les forces de frottement visqueux ne peuvent etre negligees, ou encore pour des 
ecoulements quasi permanents avec un travail utile. 

Remarques diverses : 

- pour un fluide immobile, v = et la formule de Bernoulli redonne la loi de 
l'hydrostatique. 

- la formule de Bernoulli peut etre vue comme une loi de conservation de l'energie 
mecanique : elle signifie que energie cinetique, energie potentielle de pesanteur et 
pression sont interchangeables, sans pertes. 

- la pression apparait comme une energie volumique. 
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5.3 Applications 




s 
1_ 




Pa,tm 

Patm + pgh | • s 



t 

Fig. 5.1 - Chateau d'eau connecte a un tuyau ferme 
^4 P&tm 




s 



I Pa,tm 



Patm 



5 



Fig. 5.2 - Chateau d'eau connecte a un tuyau ouvert a l'atmosphere 



5.3 Applications 



5.3.1 Problemes de vidange 



Considerons un chateau d'eau de section £ dont le fond est ouvert sur une tuyauterie 



de section S, terminee par un robinet (Fig. 5.1). Soit h la hauteur entre le robinet 
et la surface libre dans le chateau d'eau. 

Si le robinet est ferme, la pression p s'exergant sur le robinet cote fluide est p a tm+pgh 
et vaut Patm cote air. Le clapet du robinet est done soumis a une difference de pression 
nette p = pgh. C'est de cette pression que Ton parle lorsque Ton dit par exemple 
que ce robinet est a une pression de 8 bars : il s'agit aussi bien d'une hauteur egale 
a 8 bars / pg. 



Si on ouvre maintenant le robinet (Fig. 5.2) 



- le fluide dans le tuyau revient a la pression p atm 

- les vitesses aux points A et S sont reliees par v = vsS via la conservation de 
la masse. En general S ^> S et par consequent Va "C v s 

- si le tuyau est a section constante, la vitesse vaut vb dans toutes les sections du 
tube. 



5.3 Applications 
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par consequent, Bernoulli applique entre A et S indique que 

1 2 1 2 

Patm + pgk + ~pV A = p atm + -pVg 

negligeable 

l'energie cinetique du fluide a la sortie est done \pv 2 s = pgh, energie potentielle de 
pesanteur. 

Nous voyons done que la relation de Bernoulli indique qu'une hauteur est toujours 
convertible en pression ou en vitesse. C'est la notion de charge que nous definirons 
un peu plus loin. 

Notons egalement que Ton peut deduire la vitesse de vidange du chateau d'eau de 
la formule ci-dessus : 



v s = ^2g~h (5.7) 
connue sous le nom de formule de Toricelli 

a Evangeliste Toricelli (1608-1647) : mathematicien physicien italien qui fut eleve, puis 
secretaire et compagnon de Galilee pendant les 3 derniers mois de sa vie. II est l'inventeur 
de la pression atmospherique, du barometre a mercure et a donne son nom au torr unite 
de pression = 1 mm de Hg. II expliqua correctement pourquoi les variations de hauteur du 
barometre a mercure etaient 14 fois plus petites que celles du barometre a eau, en terme 
de densites. II est le premier a parler de "vide", qu'ill obtint en renversant un tube rcmpli de 
liquide dans un bain contenant le memo liquide. II a developpe egalement la notion de quantite 
de mouvemcnt. 



Une derniere remarque importante : en un point M quelconque du tuyau horizontal, 
le fluide se deplace a la vitesse v = vs, en raison de la conservation de la masse. 
La formule de Bernoulli, appliquee entre M et S indique alors que la pression en M 
vaut Patmj ce quelle que so it la longueur du tuyau! Cela signifie que le chateau 
d'eau pourrait alimenter une maison distante de plusieurs milliers de kilometres. . . 
Le paradoxe vient bien sur du modele de fluide parfait, mis en defaut sur de telles 
longueurs. Nous y reviendrons dans la section suivante. 



5.3.2 Pression dynamique. Forces sur un obstacle. 

Autre application, cette fois en ecoulement externe. Imaginons un ecoulement pa- 
rallele non perturbe a la vitesse Voo et plagons un obstacle dans cet ecoulement 



(Fig. 5.3 ). Les lignes de courant venant de 1'infini amont de l'obstacle vont le contour- 



ner de chaque cote. Intuitivement on voit qu'il existe une ligne de courant s'arretant 
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5.3 Applications 



perpendiculairement sur l'obstacle en un point A de vitesse nulle, de par la condi- 
tion d'etancheite. Un tel point est appele point d'arret, ou de stagnation de 
l'ecoulement. 

F z 




Fig. 5.3 - Ecoulement autour d'un obstacle. A est le point de stagnation ou d'arret. 
L'obstacle sub it une force de trainee F x et une force de portance F z . 

Si nous appliquons Bernoulli entre ce point et un point M m infiniment eloigne de 
l'obstacle, nous obtenons (en negligeant la pesanteur) : 



Poo + 2~P v °o 2 = Pa (5.8) 

La pression pa est appelee pression de stagnation ou pression dynamique. Elle 

est necessairement superieure a la pression infinie amont, et ce d'autant plus que la 
vitesse de l'ecoulement est elevee. On peut la mesurer avec un tube de Pitot dont 
l'element principal est un petit tube tres fin place face a l'ecoulement. Le fluide ne 
penetre pas dans le tube de telle sorte que son extremite est un point d'arret, done 



a la pression pa de la relation (5.8). 



Par ailleurs, on comprend aisement que l'obstacle subit une pression additionnelle, 
de par le mouvement du fluide, egale a l/2p\ oc 2 (e'est cette force que Ton sent 
lorsque Ton passe la main par le vitre d'une voiture). La force subie par l'obstacle 
de la part de l'ecoulement sera done de l'ordre de 



\pVoo 2 x S 

oil S est une section caracteristique de l'obstacle (par exemple 4nR 2 dans le cas d'un 
obstacle spherique). 

On decompose usuellement cette force en une composante F x parallele a l'ecou- 
lement, appelee trainee, et une force perpendiculaire Fz appelee portance (oui, 
e'est bien celle qui porte les avions . . .). Leur calcul depasse le cadre de ce cours, 
mais puisque nous avons un ordre de grandeur de ces forces, on les definit par des 
nombres adimensionnels C x et C z appeles respect ivement coefficient de trainee 
et coefficient de portance : 



5.3 Applications 
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C x 

c z 



\pVjS 



A titre indicatif le C x d'une voiture de course est de l'ordre de 0.3 alors que ce- 
lui d'une 2 CV atteint 0.6. La notation C x a donne son nom a l'une des voitures 
pionnieres en matiere d'aerodynamique dans les annees 70. . . 

Le rapport C z /C x est connu sous le nom de "finesse" en aviation et definit la capacite 
d'un avion a planer. Elle est de l'ordre de 40 pour un planeur, de 25 pour un gros 
avion de transport, et de 8 pour un petit avion de tourisme. 

Exercice : calculer l'ordre de grandeur des forces aerodynamiques sur un airbus A320 
(surface totale des ailes 122.6 m 2 , masse maximale 54 400 kg, vitesse de croisiere 
955 km/h, et masse volumique de l'air a 10000 m d'altitude = 0.41 kg /m 3 ). On 
comparera au poids (masse maximale = 54 400 kg). Meme calcul sur une voiture de 
1500 kg a 100 km/h (on intuitera les dimensions). 



5.3.3 Notion de charge 



On appelle charge ("head" en anglais) la grandeur 

H = z+>L + ll (5 . 9) 
99 2 g 

homogene a une hauteur 

— est appelee hauteur piezometrique ou charge de pression 

99 

v 2 

— est appelee hauteur capable ou charge dynamique 

2g , 

On definit aussi parfois la charge sous forme de pression : 

P = pgz + p + ]-pv 2 (5.10) 



Cette definition convertit la pression et la vitesse en une hauteur ou en une pression 
equivalentes. Ainsi, dans les hypotheses precitees : 



La formule de Bernoulli predit que la charge est constante le long d'une ligne 
de courant 
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5.3 Applications 



Cela est bien sur faux dans l'absolu a cause des frottements visqueux et l'idee d'un 
chateau d'eau alimentant une maison distante de plusieurs milliers de kilometres 
doit etre oubliee dans la pratique. 

Les differentes composantes de la charge peuvent etre representees graphiquement. 
Prenons tout d'abord le cas d'une canalisation de section constante et montante 



(Fig. 5.4). La section constante impose que la vitesse reste egalement constante le 



long du tuyau. Tragons par rapport a une reference arbitraire, les courbes devolution 
de chaque terme de la charge : z est simplement l'altitude du tuyau, z + p/pg est 
appelee ligne piezometrique, H est la ligne de charge (horizontale si la formule 
de Bernoulli est applicable), et la charge dynamique v 2 /2g est la difference entre la 
ligne de charge et la ligne piezometrique. Ici cette derniere est egalement horizontale, 
puisque la vitesse est constante. 

Ligne de charge 



H 



v 



Ligne piezometrique 



Reference arbitraire 



v 

2g 




Fig. 5.4 - Ligne de charge et ligne piezometrique pour une canalisation de section 
constante en fluide parfait. 



La figure |5.4| montre les interpretations physiques des lignes piezometriques et de 
charge. La premiere est la hauteur d'eau dans un tube connecte sur la paroi du 
tuyau. La seconde est la hauteur d'eau dans un tube relie a un embout de tube 
Pitot plonge au coeur de l'ecoulement Ce sont des consequences de la formule de 
Bernoulli (on pourra le montrer a titre d'exercice). 



Reprenons maintenant l'exemple du chateau d'eau sur le tuyau ouvert (Fig. 5.5). 
Nous prenons la reference au niveau du tuyau. La pression est constante egale a 
la pression atmospherique dans tout l'ecoulement, et la hauteur h est entierement 
convertie en charge dynamique. Notons que dans l'approximation fluide parfait, 
le fluide ne monterait pas dans un quelconque tuyau connecte sur le cote de la 



5.4 Pertes et gains de charge 
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canalisation 



Ligne de charge 



i iPatm N ^ 





\ 



\ 



2^ 



\ 



H 



h 



\ 



\ 



\ Ligne piezometrique 



m 



FlG. 5.5 - Ligne de charge et ligne piezometrique pour le probleme du chateau d'eau 
ouvert sur une canalisation a l'atmosphere. 



5.4 Pertes et gains de charge 
5.4.1 Formule de Bernoulli generalised 

La formule de Bernoulli, mal utilisee, peut fournir des resultats totalement absurdes, 
confinant parfois au comique. Tout d'abord, dans la mesure ou elle decrit un mouve- 
ment non dissipatif, il est probablement possible de lui faire inventer le mouvement 
perpetuel (rappelons le chateau d'eau alimentant une maison a quelques milliers de 
km) et de faire plaisir aux petroliers pour transporter le petrole dans les oleoducs a 
peu de frais ! Nous touchons ici aux pertes de charges occasionnees par le caractere 
visqueux du fluide, qui ne peut etre neglige dans certains cas. 

La formule de Bernoulli echoue egalement dans un autre cas fondamental : lorsque 
l'on preleve ou apporte du travail au fluide en ecoulement a l'aide d'une machine 
tournante. Dans ce cas, ou bien on ajoute de la charge a l'ecoulement, on bien on 
lui en preleve. Notons que cette configuration est loin d'etre anecdotique puisqu'elle 
est a la base des pompes et des moulins, utilises depuis l'antiquite. 

Enfin, elle est inutilisable dans certains cas un peu plus subtils, pour lesquels l'hy- 
pothese regime permanent doit etre remise en cause. 

II est done pertinent de mentionner ici les limites de la formule de Bernoulli apres 
l'avoir exposee dans les sections precedentes. Nous allons done calculer la variation 
de la charge entre la section d'entree d'un tube de courant et sa section de sortie 
dans un cas tres general. Elle est definie par : 



2 Pourquoi ? 
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5.4 Pertes et gains de charge 



, p lv 2 \ ( p lv 2 
H s -H e = ( * + — + 7T— -lz + — + -- 
99 2 g J V 99 ^9 



Pour calculer cette variation, utilisons le theoreme de l'energie cinetique (4.27) eta- 
bli au chapitre precedent, ecrit ci-dessous en regime permanent, et pour un fluide 
incompressible (divv = 0) : 



P + pgz + p—j - \ P + pgz + p—j =-^r-^r (5.H) 

ou V est le debit volumique dans le tuyau, W la puissance utile, fournie ou prelevee 
par les pales et la puissance des forces interieures visqueuses. La variation de 
charge proprement dite s'ecrit : 



+ z + —\ ~ (- + Z+ — ) ^ (5.12) 



99 ' ZgJs \P9 ZgJe pgV pgV 

On remarqueraj^] que le theoreme de l'energie cinetique fait intervenir directement 
la difference des charges entre la sortie et V entree ! On retrouve done la formule de 
Bernoulli seulement lorsque $„ = (fluide parfait) et W u = 0. 

On retiendra que : 



La difference de charge entre deux sections d'un ecoulement en tuyauterie 
s'ecrit : 

H s -H e = ^-^ = h u -h v (5.13) 
pgV pgV 

oil ho > est liee aux pertes par frottements, et h u > ou < est liee a la 
puissance mecanique utile exercee sur le fluide. 

En l'absence de machine et en fluide parfait, on retrouve la relation de Bernoulli : 

H s = H e 



Nous allons examiner individuellement les deux cas. 



3 Ce qui n'est pas etonnant puisque cette formule, comme celle de Bernoulli, decoule de la 
conservation de la quantite de mouvement. 



5.4 Pertes et gains de charge 
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5.4.2 Frottement visqueux : $ w ^ 



Nous avons deja mentionne que le terme de puissance des forces visqueuses in- 
terieures au fluide etait toujours positive. Cette energie, comme nous l'avons vu 
d'apres la formule (4.28) est soit dissipee sous forme de chaleur, soit augmente la 
temperature du fluide. Dans tous les cas, c'est une perte d'energie mecanique pour 
le fluide, et done de charge, comme en temoigne l'equation (5.12). 

Les pertes de charges sont de deux types : 



5.4.2.1 Pertes de charge regulieres 



Appelees aussi pertes de charge lineiques, elles sont liees a l'adherence du fluide sur 
les parois de la tuyauterie. Par un raisonnement dimensionnel, on montre que la 
perte de charge peut s'exprimer en fonction d'un facteur adimensionnel / appele 
coefficient de frottement : 



/ 



D2g 



(5.14) 



ou L est la longueur entre section d'entree et de sortie, et D diametre du tube. 



Le coefficient de frottement est fonction du nombre de Reynolds Re = pvD/fj,, 
du rapport L/D, et eventuellement de la rugosite du tuyau pour des ecoulements 



turbulents. En regime laminaire (Re < 2000), on montre que (Bird et al. 1960 



Chassaing, 2000, voir aussi chapitre 6 



./ 



64 
Re' 



(5.15) 



Pour des nombres de Reynolds plus eleve l'ecoulement devient turbulent. Pour cal- 
culer des ordres de grandeurs lorsque Ton a des Reynolds eleves, et un tube lisse, on 



pourra utiliser la formule de Blasius (Bird et al. , 1960) : 



/ 



0.3164 



pour 2.1 x 10 3 < Re < 10 E 



(5.16) 



Toujours pour des regimes turbulents, en tube lisse ou rugueux, on pourra utiliser 
la formule implicite 
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5.4 Pertes et gains de charge 



ou e est la taille caracteristique des asperites du tube. Cette formule presente l'in- 
convenient d'etre implicite, mais fournit des resutats continus sur l'ensemble du 
domaine turbulent. 

De fagon generale, on pourra utiliser le celebre diagramme de Moody qui fournit 
le coefficient de perte de charge lineique dans tous les cas. On retiendra que la perte 
de charge est proportionnelle a la longueur de tube considereef] 



5.4.2.2 Pertes de charge singulieres 



Elles sont liees a des "accidents" ponctuels sur la tuyauterie : elargissements ou retre- 
cissements brusques, coudes, Us, vannes, robinets, clapets, tes. Ces pertes de charge 
s'expriment egalement a partir d'un coefficient de perte de charge adimensionnel e v : 



v 2 

K = e v — (5.18) 
Le facteur e v peut etre calcule explicitement dans quelques cas peu nombreux, mais 



on a generalement recours a des formules empiriques ou a des abaques (Idelcik 



1969) 



On pourra retenir les ordres de grandeurs suivants : 

- Elargissement brusque : e v = ( 1 — j , en reference a la vitesse amont (dans 



( P8D ), 



Retrecissement brusque : e v = 0.45 ^1 — — ) , en reference a la vitesse aval, 

Tube plongeant dans un bac, e v = 1, 
Coude, e v = 0.8, 
Vanne ouverte, e v = 1.2, 



5.4.2.3 Exemple 

Nous allons reprendre l'exemple du chateau d'eau en prenant en compte les pertes 
de charge lineiques : pour cela remarquons que la perte de charge impose une de- 



croissance lineaire de la ligne de charge (Fig. 5.6). Comme altitude et vitesse sont 
constantes le long du tuyau, c'est la pression qui chute le long de l'ecoulement : ainsi, 
la ligne piezometrique reste parallele a la ligne de charge dans toute la conduite. La 
pression a l'entree du tuyau est done superieure a p atm : le ffuide est retenu par cet 



exces de pression. II apparait done clairement sur la figure (Fig. 5.6) que la charge 



dynamique, et done la vitesse, est plus faible en presence de pertes de charges. 



4 Certains goutteurs pour Farrosage sont bases sur ce principe : plusieurs dizaines de metres de 
tuyau fin sont enroules sur eux memes, de telle sorte que Feau sort a tres faible vitesse. 
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Ligne de charge en fluide parfait 




L 

Fig. 5.6 - Ligne de charge et ligne piezometrique pour le probleme du chateau d'eau 
ouvert sur une canalisation a l'atmosphere, en presence de pertes de charge. 



Cette perte de vitesse peut etre evaluee en ecrivant les variations de charge entre les 



differents points de l'ecoulement (Fig. 5.6) 



Entre A et B, il n'ya pas de pertes done Ha = Hb- Entre B et C, on a la perte de 
charge lineique done Hs — Hb = — h v soit par consequent 



Hs = Ha — h v 

soit : 



pg 2g pg 

et la vitesse d'ecoulement s'ecrit done : 



P * tm + — = ^ + h-h v (5.19) 



\llg(h - h v ) 



On constate bien qu'elle est plus faible que la vitesse predite par la formule de 



Toricelli (5.7). En comparant cette derniere a l'expression (5.7), on voit que tout se 



passe comme si le niveau dans le chateau d'eau etait diminue de la hauteur h v . On 
remarquera que si la perte de charge atteint la hauteur h du chateau d'eau, le fluide 
ne peut plus s'ecoulerj^Jla perte de charge etant proportionnelle a L, ceci se produira 
si la tuyauterie depasse une certaine longueur, ce qui montre bien 1'impossibilite pour 
un reservoir d'alimenter des maisons trop distantes. 



5 Comment y remedier? indice : lire la section suivante et/ou regarder la formule (5.121. 
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5.4 Pertes et gains de charge 



Remarque importante : il convient de bien comprendre que ce n'est pas la vitesse 
(ou le debit) qui chute le long de la conduite, mais la pression ! En revanche la 
perte de charge conduit bien a une diminution de vitesse par rapport au fluide 
parfait. 



5.4.3 Cas des machines tournantes : W u ^ 

Lorsqu'une machine contenant des pales tournantes est placee entre l'entree et la 
sortie, elle fournit de de la puissance utile au fluide (W u > 0, pompe, ventilateur) 
ou en regoit (W u < 0, eolienne, moulin a eau, turbine hydroelectrique) . Dans le 
premier cas, la charge augmente, dans le second, elle diminue. Ainsi, une pompe 
permettra par exemple de monter du fluide d'une hauteur z e a une hauteur z s en 
augmentant sa charge. A 1 'inverse, une turbine hydroelectrique preleve une partie 
de l'energie potentielle de pesanteur du fluide pour faire tourner une dynamo. La 
charge restante est essentiellement une charge dynamique (les cascades generalement 
visibles sous les barrages). 

II est important de comprendre que cette perte ou ce gain de charge est independant 
du caractere visqueux ou non du fluide. 



Un exemple d'application est presente Fig. 



5.7 



: on place une pompe {W u > 0) entre 



deux tuyauteries de meme section. La difference de charge entre entree et sortie 
s'ecrit : 



if 2 = H x + — }- = H x + h u 

pgV 

avec h u > 0. La conservation du debit implique que les vitesses sont egales en entree 
et en sortie (la charge dynamique est done constante). De plus, entree et sortie sont 
a la meme altitude et l'equation de charge s'ecrit done : 



Fa 
pg 



El 
pg 



+ K 



ce que reflete bien la difference d'elevation dans les deux tubes piezometriques. 
Conclusion : on a augmente la pression dans l'ecoulement. Cette augmentation de 
pression peut ensuite etre exploitee pour elever le fluide par exemple. 



Pour finir, les figures 5^ et |5.10| montrent des exemples de roues de moulins a 
eau, de turbines de barrages hydroelectriques et de pompes. 



Exercices de reflexion : pourquoi la formule de Bernoulli (if s = if e ) appliquee sur une 
ligne de courant entre l'entree et la sortie d'une machine n'est-elle pas applicable ? 
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Fig. 5.7 - Augmentation de la pression dans un ecoulement de meme section en 
entree et en sortie avec une pompe. 





Fig. 5.8 - Exemples de roues de moulins a eau. 



5.4.4 Applications aux reseaux de fluide 



La formule (5.12) resume clairement le probleme du dimensionnement d'un reseau 
de fluide : par exemple, pour compenser les pertes de charge entre l'entree et la 
sortie d'une conduite, il faut apporter de la puissance utile, en d'autres termes il 
faut pomper ! Dans ce cas, on compense la perte de charge par un apport d'energie 
(W u = $„). C'est par exemple le cas du transport de petrole sur des milliers de km 
dans les oleoducs. 



Si Ton veut faire monter du fluide d'une altitude z e a une altitude z s > z e dans une 
conduite de section constante (done v s = v s ), et a la meme pression, il faut ici encore 




Fig. 5.9 - Exemples de turbines de barrages. 
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5.4 Pertes et gains de charge 



fournir de la puissance utile, meme si les pertes de charge sont negligeables, done 
il faut encore pomper (W u = pgV{z s — z e ) > 0). C'est ainsi par exemple que Ton 
pompe de l'eau dans un puit. 

A l'inverse, on peut profiter de l'altitude d'un liquide (z e > z s ) pour fournir de la 
puissance utile a une turbine (W u = pgV(z e — z s ) < 0). Une surpression a l'entree 
peut egalement fournir du travail : c'est par exemple le cas de la soupape de la 
cocotte-minute qui tourne. 



A l'echelle d'une installation en genie des precedes, il est important de prendre 
en compte toutes les pertes de charge, lineiques et singulieres et de 

bien dimensionner la ou les pompes pour faire circuler le fluide avec le debit 
voulu. Lorsque la pompe est sous-dimensionnee, le fluide cesse de circuler ou 
bien revient meme de la sortie vers l'entree de la pompe. On dit alors qu'elle 
refoule. 




Fig. 5.10 - Exemples de pompes volumetriques : de haut en bas et de gauche a 
droite : pompe a engrenages exterieurs ; pompe a membranes ; pompe peristaltique ; 
pompe a palettes libres. 



Chapitre 6 
Fluides visqueux 



6.1 Introduction 



Le titre de ce chapitre est abusif dans la mesure ou tous les fluides sont en realite 
visqueux. La question est de savoir si nous devons prendre en compte ou non l'effet 
des forces visqueuses. 

Nous commencerons par expliciter la contrainte visqueuse dans le cas des fluides dits 
"newtoniens" et en remplagant cette expression dans les equations du mouvement, 
nous obtiendrons les equations dites de "Navier-Stokes". 



6.2 Le modele de fluide newtonien 



6.2.1 Approche par l'experience de Couette 



Reprenons l'experience de Couette decrite page 46 Rappelons que l'experimentation 
montre l'apparition d'un profll de vitesse lineaire au bout d'un certain temps, et que 
la force qu'il faut exercer sur la plaque pour la deplacer a une vitesse constante Uq 
est proportionnelle a la variation de la vitesse dans la hauteur du fluide : 



F v dv x 
~S = -^ ex 

Puisque la plaque bouge a vitesse constante, cette force est done opposee a celle 
qu'exerce le fluide sur la plaque. Effectuons maintenant le bilan des forces sur une 
tranche de fluide comprise entre la plaque superieure et une ordonnee z quelconque. 
On voit que : 

- suivant z, il n'y a pas de mouvement, les seules forces sont les forces de pression, 
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6.2 Le modele de fluide newtonien 



nous allons retrouver la loi de l'hydrostatique, comme si le fluide etait immobile. 
- 1' acceleration du fluide etant nulle (nous y reviendrons) , la force tangentielle exer- 
cee par la plaque sur la tranche de fluide est egale a celle qu'exerce cette tranche de 
fluide sur la couche de fluide inferieure. Par consequent, la contrainte tangentielle 
exercee en tout point par le fluide superieur sur le fluide inferieur est constante 
dans l'epaisseur et vaut : 

dv x 

La contrainte normale exercee par la couche superieure sur la couche inferieure est 
celle de pression et nous pouvons done ecrire la contrainte totale sous la forme : 



Remarquons que le vecteur normal sortant a la couche de fluide inferieure est e z , 
et la relation ci-dessus nous fournit la troisieme colonne du tenseur des contraintes 
(qui rappelons-le, est la force exercee par unite de surface normale a Oz). De plus 
nous savons que la force de pression est isotrope (elle est toujours orthogonale a la 
surface considered) . Le tenseur des contraintes s'ecrit done : 



-P 



1 1 



dv x 
dz 



-p 



? _ 



V 



Dans la direction y, il n'existe pas d'autre force que celle de pression (par symetrie 
de l'ecoulement), par consequent la deuxieme colonne est nulle en dehors du terme 
—p. Enfin examinons la force tangentielle exercee sur une surface perpendiculaire 
a l'ecoulement : l'existence d'une composante suivant y imposerait un mouvement 
dans cette direction, et n'existe done pas, par consequent a xy = 0. Enfin, on peut 
demontrer tres generalement que le tenseur des contraintes est forcement symetrique, 
et nous obtenons finalement : 



-p 



/I 



dv x 
dz 



-p 

dv x 



(6.1) 
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6.2.2 Equations du modele newtonien 

L'experience de Couette nous montre que la force visqueuse nait de l'existence de 
couches de fluide voisines se deplagant a des vitesses differentes. L'expression ^ 
est une derivee spatiale d'une composante de la vitesse. Voyons combien de derivees 
spatiales de la vitesse nous pouvons definir dans le cas general : il y a trois com- 
posantes u, v, w que nous pouvons chacune deriver par rapport aux trois variables 
d'espace x,y,z, ce qui fait 9 derivees. On les regroupe dans le tenseur gradient, 
qui s'ecrit, en coordonnees cartesiennes : 



gradv 



dv x 


dv x 


dv x 


dx 


dy 


dz 


OVy 


dv y 


dVy 


dx 


dy 


dz 


dv z 


dv z 


dv z 


dx 


dy 


dz 



(6.2) 



II est vraisemblable que le tenseur des contraintes visqueuses pour un fluide newto- 



nien fasse intervenir ce tenseur. En examinant le tenseur (6.1) obtenu pour l'expe- 



rience de Couette, il est facile de voir que la relation est la suivante : 



cr = — pi + fx I T gradv + gradv J (6-3) 



Un fluide dont suivant cette loi sera dit newtonien. Attention, beaucoup de fluides 
ne sont pas newtoniens. Citons en vrac les pates et suspension de fagon generale, la 
plupart des fluides alimentaires (creme, yaourt, ketchup), le sang, les polymeres a 
l'etat liquide . . . 



6.3 Equations de Navier-Stokes 



II nous est maintenant possible d'expliciter le terme visqueux dans toutes les equa- 
tions presentees jusqu'alors. Nous nous cantonnerons au fluide incompressible. Ainsi 
l'equation locale (4.31) devient, en ecrivant div<r„ a l'aide des formules d'analyse 
vectorielle rappelees en annexe [B] : 



<9v 

P-qI + p(v.V)v 



gradp + pg + /A7 2 v 



(6.4) 
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6.3 Equations de Navier-Stokes 



Cette equation, associee a la conservation de la masse 

divv = 0, (6.5) 
constitue les equations de Navier-Stokes incompressible]]] 

6.3.1 Adimensionnalisation 

Dans tout ecoulement, on peut trouver une longueur de reference : le diametre d'un 
tuyau, la longueur d'une plaque, le diametre d'une sphere, la dimension transverse 
moyenne (perpendiculaire a l'ecoulement) d'un vehicule. De meme, on peut trouver 
une vitesse de reference : vitesse incidente sur un obstacle, vitesse calculee a partir 
d'un debit. . .. Posons L cette longueur et V cette vitesse, et rapportons toutes les 
variables spatiales ainsi que le champ de vitesse respectivement a la longueur et la 
vitesse de reference (on notera les variables adimensionnelles avec une etoile) : 

x = Lx* y = Ly* z = Lz* 
v x = Vv* v y = Vv* v z = Vv* z 

On rapporte egalement le temps a un temps de reference T : ce peut etre par exemple 
une periode pour un ecoulement oscillant ou un temps caracteristique de l'application 
d'une contrainte externe a un ffuide. On pose 



t = Tt* 

Enfin, on peut, pour des applications bien specifiques, ramener la pression a une 
grandeur de reference po, par exemple la pression atmospherique 



V = PoP , 

mais dans la plupart des cas, on prendra comme pression de reference po la pres- 
sion dynamique (voir le chapitre sur Bernoulli) pV 2 . En remplagant les variables 



originales par ces variables dans l'equation de Navier-Stokes (6.4), on obtient 



1 Qui constitue le modele theorique de base pour la plupart des ecoulements. La plupart des 
ecoulements compressibles peuvent etre traites par les equations d'Euler, en negligeant done le 
caractere visqueux du fluide. Les equations de NS compressibles sont utilisees dans un nombre de 
cas assez reduit. Notons par aillcurs que rexistence et l'unicite de solutions aux equations de NS 
constituent l'objet de l'un des dix problemes proposes par la fondation Clay, recompense par un 
million de dollars. . . ce qui n'empeche pas le chercheur et l'ingenieur d'utiliser ces equations au 
quotidien dans un large cvcntail d' applications. II est interessant de noter qu'il s'agit pourtant d'un 
problcmc de mccaniquc classique, qui n'a pas encore livre tous ses mysteres. 
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Vdv* V 2 . , po gradV ^„*2 * 

I v .V v = - — — + £ H V v , 

T dt* L pL p L 2 



soit, en multipliant toute l'equation par L/V 2 , 



St ^ + (v*.V*)v* = gradV + Fr^ + -*-V*V, 
at* Eu j Re 



(6.6) 



avec les nombres adimensionnels : 
LV 



Reynolds 


Re 


Strouhal 


St 


Euler 


Eu 


Froude 


Fr 



Inertie du fluide 



v Forces visqueuses 
L Temps caracteristique du mouvement du fluide 

VT Temps caracteristique de la sollicitation 

pV 2 Pression dynamique 

Po Pression statique 

gL Poids 

V 2 Inertie 

Ces nombres adimensionnels sont tres importants car leur calcul fournit deja des 
informations tres importantes sur l'ecoulement, sans avoir a resoudre les equations de 
NS. Le plus important de tous est le nombre de Reynolds qui quantifie l'importance 
de l'inertie du fluide par rapport aux forces visqueuses. 



6.3.2 Classification des ecoulements 

Les ecoulements se comportent de fagon tres differente selon la valeur du nombre de 
Reynolds. 

- Pour des nombres de Reynolds tres faibles, l'ecoulement est domine par les forces 
visqueuses. Visuellement, on observe un ecoulement tres "penible", que Ton ima- 
gine difficile a mettre en mouvement, et facile a arreter. Ce type d'ecoulement est le 
moins frequent dans la nature, mais intervient par exemple dans des problemes de 
lubrification, dans des suspensions de petites particules. On parle d'ecoulement 
de Stokes ou ecoulement rampant. Nous verrons que ce type d'ecoulement est 
prevu par la reduction des equations de Navier-Stokes dans la limite Re <^ 1. 

- pour des Reynolds de l'ordre de l'unite, l'ecoulement se "complique", par exemples 
des recirculations apparaissent en aval des obstacles ou dans les singularites d'un 
tuyauterie (coudes, retrecissements, . . .). L'ecoulement reste malgre tout station- 
naire et laminaire sauf bien sur s'il est sollicite de fagon instationnaire. Pour des 
ecoulements autour d'obstacles symetriques, on note une dissymetrisation amont- 
aval. 

- Pour des Reynolds encore superieurs, des variations temporelles periodiques ap- 
paraissent dans le champ de vitesse. 
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6.4 Ecoulements visqueux classiques 



- Pour des Reynolds tres grands, l'ecoulement est completement destabilise , soit 
dans une certaine zone de l'espace (dans la couche limite et le sillage des obstacles), 
soit dans tout l'espace (dans un tube) et au champ de vitesses principal (obser- 
vable macroscopiquement) se superposent des fluctuations aleatoires a plusieurs 
echelles : c'est la turbulence. 

La classification precedente est tres reductrice, revolution de l'ecoulement dependant 
de la geometrie et des conditions aux limites. Cependant on retrouve des points 
communs dans une meme famille d'ecoulements et quelques cas limites importants 
peuvent etre traites de fagon analytique : 

- les ecoulements rampants (a faible Reynolds) : ils permettent le calcul des forces hy- 
drodynamiques dans de nombreux cas, et debouchent sur la theorie des suspensions, 
des milieux poreux, ... Les equations de Navier-Stokes se simplifient grandement car 
le terme (v.V)v se simplifie. 

- les problemes de couches limites laminaires, utilisables tant que l'ecoulement n'est 
ni destabilise, ni decolle de la paroi (comme c'est le cas pour les profils epais style 
cylindres ou spheres). On peut calculer rigoureusement le champ de vitesses dans la 
couche limite. 

- En regime turbulent, les problemes de couche limite en ecoulement ouvert, et de 



couche turbulente a la paroi des tuyaux admettent des solutions approchees. Guyon 



et al. (2001) 



Nous traiterons done ces problemes separement. Notons tout d'abord qu'un certain 
nombre de problemes "classiques" peuvent se resoudre analytiquement, en raison 
de la simplicity de la geometrie. Nous allons aborder certains de ces problemes des 
maintenant. 



6.4 Ecoulements visqueux classiques 



6.4.1 Ecoulement de Couette 



Nous sommes desormais familiers avec cet ecoulement. Essayons de retrouver les 
resultats precedents d'apres les equations de Navier-Stokes. 

Notons tout d'abord que le champ des vitesses est oriente suivant x. On a done 
v = [v x {x, y, z),0, 0]. L'equation de conservation de la masse fournit alors 



JT = ° 
ox 

Done v x ne depend que de y et z. On suppose l'ecoulement d'extension infinie suivant 
y et par consequent v x est independant de y. Le champ de vitesse s'ecrit done 
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finalement v = (v x (z), 0, 0). 

Remarquons alors que le terme (v.V)v est identiquement nul. En effet le produit 
scalaire formel (v.V) se reduit a v x J^ et la vitesse ne depend pas de x. Ce terme 
disparait done a priori tant que l'ecoulement reste dans la direction x, sans hypothese 
sur le nombre de Reynolds. L'ecoulement etant stationnaire, l'equation de Navier- 
Stokes se reduit done a : 



= — gradp + fiV 2 v + pg (6.7) 
soit en projection sur les 3 axes : 



dp | <Pv x 

dx dz 2 
dp 



dy 
dp 

dz 



(6.9) 

99 (6.10) 



ou nous avons ecrit un d droit dans la premiere equation, puisque v x ne depend que 
de z. La deuxieme equation nous indique que p ne depend que de x et z. La troisieme 
equation s'integre simplement en 



p(x, z) = —pgz + K{x) 

qui indique que la pression varie comme en hydrostatique suivant la direction de la 
pesanteur, independamment de l'ecoulement. Par ailleurs, aucun gradient de pression 
parallelement aux plaques n'est applique a l'ecoulement, et la pression est done 
independante de x, [j done K(x) = K. La premiere equation s'ecrit done 

d 2 v x 

^ = (6.11) 

et on obtient done, avec les conditions aux limites 1^(0) = et v x (h) = U , un profll 
de vitesses lineaire : 



v x (z) = U ~ (6.12) 
h 

2 Les ecoulements de Couette reels sont realises dans des canaux annulaires. Dans ce cas on 
concoit aisement que la pression ne puisse pas varier le long de l'ecoulement, puisque ce dernier se 
refcrme sur lui-mcme. . . 
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6.4 Ecoulements visqueux classiques 



6.4.2 Ecoulement de Poiseuille 
6.4.2.1 Calcul de l'ecoulement 



Dans cet ecoulement, le moteur de l'ecoulement n'est plus un cisaillement (done 
la force visqueuse), mais un gradient de pression entre l'entree et la sortie de la 
conduite. On peut par exemple imposer une surpression a l'entree avec une pompe 
ou un chateau d'eau, ou une depression a la sortie par une aspiration. Dans tous 
les cas on notera p e et p s les pressions respectivement a l'entree et a la sortie de la 
conduite de longueur L. 

Nous traiterons directement le cas d'un tuyau cylindrique d'axe z. On note r et 9 
les directions radiales et orthoradiales. Le champ de vitesse est dirigee uniquement 
suivant z done v g = et v r = et done v = [0, 0, v z (r, 9, z)]. 

L'equation de conservation de la masse divv fournit alors 



On suppose de plus que l'ecoulement est symetrique par rapport a l'axe du tube 
et done v z ne depend finalement que de r. II est facile de voir que comme pour 
l'ecoulement de Couette, le terme (v.V)v est identiquement nul, et en negligeant de 
plus la pesanteur, les equations de Navier-Stokes se reduisent a 









dp 

dr 
dp 

dp Id 

dz ^ r dr 



dv z 
dr 



(6.13) 
(6.14) 
(6.15) 



La pression ne depend done que de z, et la troisieme equation s'ecrit 



dp 

dz 



1 d 



dv 7 



r dr \ dr 



Dans cette equation le membre de gauche est une fonction de z seulement, et celui 
de droite est une fonction de r seulement. Leur egalite Vr, Vz indique done que les 
deux termes sont independants de r et de z, et sont done egaux a une constante A. 



L'equation (6.11) revient done a un systeme de deux equations differentielles : 
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dp 

dz 



1 d ( dv 



r dr \ dr 



A, 
A 



(6.16) 
(6.17) 



La pression varie done lineairement le long de l'ecoulement, et avec les conditions 
aux limites en pression a l'entree et a la sortie, on obtient la constante A : 



A 



Ps-Pe 



(6.18) 



Une premiere integration de l'equation (6.17) fournit alors 



dv z _ A B 
dr 2/j, r 



(6.19) 



la constante d'integration B est forcement nulle sans quoi la vitesse deviendrait 
infinie sur l'axe du tube. En utilisant ensuite la condition d'adherence sur la paroi 



du tube v z (R) = 0, l'integration de (6.19) fournit : 



v z (r) = A p _ ^) 



soit en explicitant la constante A grace a (6.18) : 



v z (r) 



Pe-Ps 

4/iL 



(R 2 - r 2 ) 



(6.20) 



On voit done que le profil de vitesses est parabolique, maximum sur l'axe, et que 
v z > si p e > p s . L'ecoulement est done oriente des fortes pressions vers les faibles. 



6.4.2.2 Calcul du debit. Perte de charge en ligne 

On peut en deduire le debit a travers le tube en integrant le profil de vitesse sur la 
section droite : 



Q — I v z (r)2iTr dr = — — —2tt / (R 2 — r 2 )r dr 

4/iL J Q 



L'integrale se calcule facilement avec le changement de variable y = (R 2 — r 2 )/R 2 
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J\r 2 - r 2 )r dr = J° R 2 y {^\ R ^j d V 



4 



On obtient done finalement le debit, en fonction du rayon ou de la section du tube : 

= ( Pe - Ps )nR* = (p e -p s )S 2 
^ 8/jL 8nfiL v ; 

On peut en deduire la vitesse moyenne : 

Pour une vitesse moyenne donnee v m , on peut done en deduire la chute de pression 
le long de l'ecoulement, qui est aussi la perte de charge en ligne pour un tuyau 
horizontal de section constante : 



8/j,Lv m 

Pe~Ps = 

En ramenant cette chute de pression a la pression dynamique 1/2/w^, et en intro- 
duisant le nombre de Reynolds Re = pv m D / 1 \i on obtient : 



Pe -p s _ 16 f±L _ 64 L_ 
\pv 2 m ' pR 2 v m ReD 



On voit done, comme annonce dans la section [5. 4. 2[ que la perte de charge en ligne 
est donnee par : 

Ap = f 2 PVm D 

avec / coefficient de frottement = 64/Re. Ce dernier resultat est valable tant que 
rhypothese de mouvement seulement suivant z est verifiee, autrement dit que le 
mouvement est laminaire. Dans la pratique, cela est verifie pour Re < 2000, valeur 
au-dela de laquelle des fluctuations turbulentes apparaissent. 



6.4.2.3 Application aux milieux poreux 



La seconde expression du debit dans (6.21) est interessante et permet de repondre 
a la question suivante : pour un meme gradient de pression impose, vaut-il mieux 



6.4 Ecoulements visqueux classiques 
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utiliser une conduite de section S ou 2 conduites de section S/2 (de fagon a garder 
une section totale constante) ? La reponse est dans l'exposant 2 de la surface S : 



dans le premier cas on trouve un debit total Qi donne par (6.21 ), dans le second un 
debit total Q 2 = <3i/4+Qi/4 = Qi/2, done deux fois inferieur. Si on prend n tuyaux 
de section S/n, on obtiendra un debit Q n = Qi/n 2 + Qi/n 2 ■ ■ ■ + Qi/n 2 = Qi/n. II 
est done plus difficile de faire couler un ffuide dans plusieurs petits tuyaux que dans 
un gros, a section totale constante. 

Calculons la vitesse moyenne v m = Q n /S pour n tuyaux de section S n = S/n a 



partir de (6.21 ) 



Qn _ Ql_ _ (ge ~ Ps) S _ S n /8n (p s -p e 

S nS 8nuL n fi L 



(6.23) 



La derniere expression est a la base de la loi de Darcy, qui decrit de fagon sim- 
plifiee 1'ecoulement dans le milieux poreux, qui a l'echelle microscopique peut etre 
vu comme un empilement de petits tuyaux de sections individuelles S n . Dans ces 
conditions (p s — p e )/L est le gradient de pression local, et la loi de Darcy s'ecrit : 



k A 

v = grad p 



(6.24) 



En comparant cette expression a (6.23), on voit que k correspond a S n /8ir et est 
homogene a une surface. Cette grandeur s'appelle la permeabilite du milieu poreux. 



6.4.3 Ecoulements unidirectionnels. 



Prenons le cas general d'un ecoulement unidirectionnel, appelons x l'axe de 1'ecou- 
lement, y l'axe perpendiculaire au plan forme par x et la pesanteur et z de fagon 
a former un triedre direct. II est facile de voir que, comme pour les ecoulements 
de Couette et de Poiseuille, le terme en (v.V)v de l'equation de Navier-Stokes est 
identiquement nul, et que la projection de cette derniere sur les 3 axes s"ecrit 









dp 

dr 
dp 

dp Id 

dz ^ r dr 



dv z 
dr 



(6.25) 
(6.26) 
(6.27) 



Chapitre 7 



Ecoulements rampants 



7.1 Equations 



En prenant la limite Re <C 1 dans l'equation (6.6), le terme (v*.V*)v* devient 
negligeable, et on obtient : 

RfiSt^ = -5^gradV + Refr- + V*V, (7.1) 
at* Eu g 

Dans la plupart des problemes, l'ecoulement est permanent et le terme instationnaire 
disparait. Si l'ecoulement n'est pas permanent, par exemple si une sollicitation exte- 
rieure est periodique (plaque vibrante ou glissante periodiquement), on peut quand 
meme negliger le terme instationnaire a condition que le produit Re St soit tres 
petit devant 1, soit St <C 1/Re, ou bien T ^> L 2 /u qui donne une borne inferieure 
du temps caracteristique T de la sollicitation. 

De fagon similaire, la pesanteur sera negligee ou non selon la valeur du produit Re Fr. 
Nous ne nous interesserons dans ce qui suit seulement a des problemes stationnaires 
ou consideres comme tels, sans influence de la pesanteur. Sous forme dimensionnelle, 
l'equation de Navier-Stokes se simplifie alors en : 



gradp + /iV 2 v = 



(7.2) 



Cette equation, comme sa forme generate (7.1) est beaucoup plus simple que l'equa- 



tion originale, car le terme non-lineaire (v.V)v a disparu. Les equations aux derivees 
partielles lineaires possedent de nombreuses proprietes mathematiques importantes, 
et un grand eventail d'outil mathematique est utilisable pour leur resolution. 
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Chapitre 8 



Couche limite 



On s'interesse dans cette partie a des ecoulements tels que S > 1, mais suffisam- 
ment petits pour que l'ecoulement reste laminaire et stationnaire. Dans une telle 
configuration, l'ecoulement a une direction privilegiee. 
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Annexe A 



Puissance du poids et premier 
principe 



II s'agit d'e valuer l'integrale : 

/ = / / / pg.v dV 




v 



sur le tube de courant. Remarquons que la pesanteur g derive d'un gradient. Si l'on 
oriente l'axe Oz vers le haut, nous avons g = grad(— gz). En utilisant de plus la 
formule div(au) = adivu + (grada).u, nous avons : 



pg.v = — div (pgzv) + gz div (pv) 

Utilisons de plus l'equation de conservation de la masse pour evaluer div (pv), nous 
obtenons : 



pg.v = - div (pgzv) - gz^ 

= -div 

en remarquant que gz n'est pas une fonction du temps. Par consequent l'integrale / 
devient : 



V V 
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La premiere integrate peut etre transformed en integrate de surface par la formule 
d'Ostrogradski, soit : 



J J pgz v.n 



dV 




d(pgz) 



dV 



pgz v.n dV 



dt 




Of 
pgz dV 



L' integrate de surface est nulle sur toutes les parois solide fixes en raison de la 
condition d'etancheite, et elle se reduit aux surfaces d'entree S e , de sortie S s et des 
pales S u . En supposant que sur les surfaces d'entree et sortie, z, p et v sont a peu 
pres constants, on obtient : 



/ = p e gz e v e S e - p s gz s v s S s - J J pgz v.n dV - — III pgz dV 

S u 




En introduisant le debit massique M = pvS, on obtient 



/ = M e gz e - M s gz s 



tit dE p 

pgz v.n dv — 

dt 



b u 



L' integrate sur S u est une contribution supplement aire de la puissan ce de s pales 
et peut done etre integree a la puissance utile W u . Ainsi, la formule (4.25) s'ecrit 
finalement 



j t (U + K + E p ) = M e (h e + V l + g Z \ -M s (h s + V l+ gZs ) +W U + Q (A.l) 



Annexe B 

Quelques formules d'analyse 
vectorielle 



v 2 



(v.V)v = (grad v).v = grad — + rot v x v (B.l) 



B.l Formules de Green et Ostrogradski 




gradadV = JJ an dS (B.2) 
\ s 

divudV = J J u.ndS (B.3) 
\ s 

div T dV = J J T.n dS (B.4) 

r s 





Consequence (obtenue a partir de ( |B.2 ) avec a = 1) : 

n dS = (B.5) 
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B.2 Derivation de produits 



B.2 Derivation de produits 



div (cm) = a div u + grad a.u (B.6) 
grad (ab) = a grad b + b grad a (B.7) 



div(T.u) = (div T).u + T : grad u (B.8) 

(le produit tensoriel A : B est le scalaire ^ijdijbij. C'est une generalisation aux 
tenseurs du produit scalaire de deux vecteurs). 



B.3 Derivee particulaire 

La variation d'une grandeur locale g(x,y,z,t) associee a un petit volume de fluide 
se deplagant d'un point de coordonnees r = (x,y,z) a un point r + dr = (x + dx,y + 
dy, z + dz) pendant un intervalle de temps dt est : 

dg = ^-dt + ^-dx + T^-dy + ^-dz (B.9) 
at ox dy oz 

Mais on a dr = \dt et par consequent : 

|^| + grad9 .v (B,0) 

On voit que la variation de g est liee a deux causes : sa variation temporelle in- 
trinseque dg/dt (par exemple l'echauffement lie a une compression locale) et sa 
variation a cause du deplacement de la particule fluide dans un champ inhomogene 
(la particule chauffe parce qu'elle voyage d'une zone froide vers une zone chaude par 
exemple). La grandeur ci-dessus est appelee derivee particulaire. On la note souvent 
Dg/ Dt pour bien specifier que g est une fonction de plusieurs variables, et ne pas la 
confondre avec la derivee classique d'une fonction d'une variable. 



Si on applique l'equation (B.10) aux 3 composantes d'un vecteur w il est facile de 



voir que la relation se generalise a l'aide du tenseur gradient en : 

dw dw 



dt 9t +s radw - v ( Rn ) 



B.3 Derivee particulaire 
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L'equation de conservation de la quantite de mouvement utilise une telle relation ou 
le vecteur w n'est autre que la vitesse v. 

Si g est une grandeur extensive massique, on demontre facilement l'identite suivante, 
en utilisant l'equation de conservation de la masse : 



Cette formule est utile pour exprimer les equations sous forme dite "conservative", 
en particulier pour la conservation de l'energie. On reconnait, dans l'interieur de 
l'operateur divergence, l'expression d'un flux convectif. 



Annexe C 

Rappel sur les forces d'inertie 

C.l Rappel : composition des accelerations 




On considere un repere TV lie a O' en mouvement par rapport a un referentiel 
galileen 1Z lie a O. On veut ecrire une relation entre les accelerations dans ces deux 
referentiels. On montre que pour un point A en mouvement quelconque : 

= a-A/iz' + a e + a c (C.l) 
ou a e est l'acceleration d'entrainement : 

a e = ao 7 w + A O'A + u n , /n A {u w in A O'A). (C.2) 

Le premier terme est lie a la translation, le second a l'acceleration angulaire, le 
troisieme a la rotation, de TV par rapport a 1Z. L'acceleration de Coriolis a c s'ecrit : 

a c = 2u w /n A VA/W- (C3) 

Remarque importante : pour un point A immobile dans le referentiel non galileen 
1Z ; , l'acceleration de Coriolis est nulle. Ce sera le cas d'un fluide immobile dans 71'. 
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C.2 Forces d'inertie 
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C.2 Forces d'inertie 



On cherche a voir dans quelle mesure on peut ecrire la loi de la dynamique dans un 
referentiel non galileen. Commengons par l'ecrire dans le referentiel galileen 1Z pour 
un point A de masse m : 



SF 



(C.4) 



soit en vertu de (C.l) 



ma.A/ii' — SF — ma e — ma c (C.5) 

On voit done que la loi de la dynamique est transposable dans TV sous reserve de 
prendre en compte deux forces supplementaires fj e et fj c dites «d'inertie d'entraine- 
ment» et «d'inertie de Coriolis». La force d'entrainement est couramment ressentie 
dans une voiture qui freine ou qui prend un virage, dans un manege forain, dans un 
ascenseur au demarrage ou a l'arret. Elle est egalement mise a profit dans les vols 
paraboliques, ou elle compense exactement la pesanteur a l'interieur de la cabine de 
l'avion. 

La force d'inertie de Coriolis est responsable notamment de la deviation vers Test 
(pendule de Foucault) et du mouvement des masses nuageuses. En revanche, contrai- 
rement a une idee regue elle n'est pas responsable du tourbillon de vidange dans un 
lavabo, et dont le sens de rotation serait pretendument oppose dans les deux hemi- 
spheres. 



Annexe D 



Demonstration du theoreme de 
l'energie cinetique 



Reprenons l'equation locale de conservation de la quantite de mouvement (4.31) en 
utilisant la formule (B.l) et faisons-en le produit scalaire par le vecteur vitesse v. 
Nous obtenons (en remarquant que le produit mixte (rot v x v).v est nul : 



dv v 2 = 

p-T^.v + pgrad —.\ = - gradp.v + pg.v + (diver „).v (D.l) 

Examinons tout d'abord le premier membre. Nous avons tout d'abord : 



dv 



d v 2 



0_ 
dt 



v 2 dp 
~2di 



Nous avons egalement, en utilisant (B.6) 



pgrad — .v = div ( p — v 
2 V — 



div (pv) 



et le premier membre s'ecrit done : 



d i v \ ,. / v 
p— + div p— v 



| + div(pv))^ 



= d'apres 
l'equation de 
conservati on de 
la masse (4.29) 
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Examinons maintenant le second membre de (D.l ). On a en utilisant (B.6) et (B.8) : 



— grad p.v+pg.v+(div cr v ).v = — div (pv)+pg.v+div (cr v .v)+p diw— cr v : grad v 



L'equation (D.l) s'ecrit done finalement : 



d ( v \ I v 

p— + div p— v 



Of 



div (pv) + pg.v + div (cr„.v) + p div v — <r v : grad v 

(D.2) 



En integrant c ette eq uatio n sur tout le volume V, en utilisant les formules de Green- 
Ostrogradski (B.3), (B.4), et en rappelant que le tenseur er est symetrique, on ob- 
tient : 




v 



d ( v 



Of 



p—\ dV+ / / p— (v.n) dS 



-pv.n dS + J J J pg.v dv + JJ (^- n )- v dS 

v s 



Puissance des 
forces exterieures 

W 




- I I I p div v dV — J J I cr v : grad v dV 

v' 




Puissance des Puissance des 

forces interieures forces visqueuses 
de pression interieures 

= -*„ 



(D.3) 



qui est bien l'equation (4.27). Nous montrerons dans le cours sur les fluides visqueux 



que le produit a v : grad est toujours positif, si bien que la puissance des forces 
visqueuses interieure est toujours negative. C'est done un puit d'energie cinetique 
pour le fluide. 



Annexe E 



Diverses formes de Pequation de 
conservation de Penergie 



En soustrayant le theoreme de l'energie cinetique (D.2) du premier principe sous 



forme locale (4.32), on obtient 



d(pu) 
dt 



+ div (pu) = —p div v + <r v : grad v — div q 



(e.i; 



Les deux premiers termes du membre de droite representent la puissance des forces 
interieures, respectivement de pression (reversible) et visqueuses (irreversible). Comme 
indique dans l'annexe precedente, la puissance des forces visqueuses interieures 
cf v : grad v est toujours positive, et nous le noterons desormais <p v . Ce terme 
contribue toujours a une augmentation de l'energie interne du ffuide, ce qui se tra- 
duit par exemple de rechauffement des fluides de lubrification, cisaille entre deux 
pieces mobiles. 

Le membre de gauche peut etre egalement ecrit pDu/Dt, ou D /Dt = d/dt + v. grad 
est la derivee particulaire, d'ou : 




ou encore, sous forme conservative, en utilisant (B.12) : 



d (pu) 
dt 



+ div (puv) = —p div v + <p v — div q 



(E.3) 



Les equations qui suivent sont toutes issues de cette derniere, mais utilisent des 
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fonctions d'etat autres que l'energie interne. On peut par exemple utiliser l'enthalpie, 
donnee sous forme differentielle par : 



h = u + 



V 
P 



(E.4) 



et par consequent 



Du _ Dh I Dp p Dp 
~D~i~ ~D~i~ ~p~D~t + p^lDt 



(E.5) 



La conservation de la masse s'ecrit par ailleurs 



Dp 
Dt 



+ p div v = 



et l'equation (E.5) s'ecrit done egalement 



Du Dh 1 Dp p 

= div v 

Dt Dt p Dt p 



(E.6) 



En remplagant cette expression dans (E.2), on obtient done 



Dh Dp , 
P Di=Di + 4>v - diyq 



(E.7) 



ou encore, sous forme conservative, en utilisant (B.12) : 



d(ph) ,. , , . Dp 

+ div (phv) = + <j> v - div q. 



(E.8) 



II convient de noter ici une simplification importante de cette equation, avec les 
hypotheses suivantes : 

- pour un gaz parfait ou un liquide incompressible dh = c p dT, 

- le flux de chaleur conductif est donne par la loi de Fourier q = — AgradT. 



L'equation (E.8) s'ecrit alors 



(E.9) 
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L'expression a l'interieur de l'operateur div contient deux flux : le flux convectif 
pCpTv et le flux conductif — AgradT. Le terme en DP/Dt est negligeable lorsque 
les variations spatiales et temporelles de la pression sont faibles, c'est souvent le cas 
pour des liquides. La puissance des forces interieures <f> v est souvent negligeable pour 
les gaz, et egalement dans de nombreux cas pour les liquides. C'est en revanche ce 
terme qui produit un echauffement du fluide dans les ecoulements de lubrification. 

La conductivity thermique peut etre sortie de l'operateur divergence si les variations 
de temperature au sein de l'ecoulement sont faibles. En effet, en general A varie 
avec T, et il augmente notablement pour les gaz chauds. 

On notera enfin que si la vitesse est nulle (fluide au repos ou solide), v, (f> v et DP/Dt 
sont identiquement nul, et on retrouve l'equation de la chaleur : 

dT 

pc p — = AV 2 T. (E.IO) 
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